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Abstract 

We consider an elliptic self-adjoint first order pseudodifferential opera- 
tor acting on columns of m complex-valued half-densities over a connected 
compact n-dimensional manifold without boundary. The eigenvalues of the 
principal symbol are assumed to be simple but no assumptions are made 
on their sign, so the operator is not necessarily semi-bounded. We study 
the spectral function, i.e. the sum of squares of Euclidean norms of eigen- 
functions evaluated at a given point of the manifold, with summation carried 
out over all eigenvalues between zero and a positive A. We derive a two-term 
asymptotic formula for the spectral function as A tends to plus infinity. We 
then restrict our study to the case when m = 2, n = 3, the operator is 
differential and has trace-free principal symbol, and address the question: 
is our operator a massless Dirac operator? We prove that it is a massless 
Dirac operator if and only if the following two conditions are satisfied at 
every point of the manifold: a) the subprincipal symbol is proportional to 
the identity matrix and b) the second asymptotic coefficient of the spectral 
function is zero. 
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1 Main results 

The aim of the paper is to extend the classical results of [TT] to systems. We are 
motivated by the following two observations. 

• To our knowledge, all previous publications on systems give formulae for the 
second asymptotic coefficient that are either incorrect or incomplete (i.e. an 
algorithm for the calculation of the second asymptotic coefficient rather than 
an actual formula). The appropriate bibliographic review is presented in 
Section [T51 

• Systems are fundamentally different from scalar operators in that spectral 
analysis of systems reveals a very rich geometric structure. An important ex- 
ample of an elliptic system is the massless Dirac operator which is examined 
in detail in our paper. 
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Consider a first order classical pseudodifferential operator A acting on columns 
v = (i>i . . . v m J of complex- valued half-densities over a connected compact n- 
dimensional manifold M. Throughout this paper we assume that m > 2 and 
n > 2. 

We assume the coefficients of the operator A to be infinitely smooth. We also 
assume that the operator A is formally self-adjoint (symmetric): J M w*Avdx = 
J M (Aw)*v dx for all infinitely smooth v, w : M — > C m . Here and further on the 
superscript * in matrices, rows and columns indicates Hermitian conjugation in 
C m and dx :— dx 1 . . . dx n , where x — (x , . . . , x n ) are local coordinates on M . 

Let A\(xX) be the principal symbol of the operator A. Here £ = (£1,... , £ n ) 
is the variable dual to the position variable x; in physics literature the £ would 
be referred to as momentum. Our principal symbol A\ is an m x m Hermitian 
matrix-function on T'M :— T*M \ {£ = 0}, i.e. on the cotangent bundle with the 
zero section removed. 

Let h^\x,$,) be the eigenvalues of the principal symbol. We assume these 
eigenvalues to be nonzero (this is a version of the ellipticity condition) but do not 
make any assumptions on their sign. We also assume that the eigenvalues (x, £) 
are simple for all (x, £ ) G T'Al. The techniques developed in our paper do not work 
in the case when eigenvalues of the principal symbol have variable multiplicity, 
though they could probably be adapted to the case of constant multiplicity different 
from multiplicity f . The use of the letter "ft." for an eigenvalue of the principal 
symbol is motivated by the fact that later it will take on the role of a Hamiltonian, 
see formula (|1.13l) . 

We enumerate the eigenvalues of the principal symbol h^(x, £) in increasing 
order, using a positive index j = 1, . . . ,m + for positive h,V'(x,£) and a negative 
index j = —1, . . . , —mT for negative ft,^' (#,£). Here m + is the number of positive 
eigenvalues of the principal symbol and m~ is the number of negative ones. Of 
course, m + + m~ — m. 

Under the above assumptions A is a self-adjoint operator, in the full functional 
analytic sense, in the Hilbert space L 2 (M: C™ 1 ) (Hilbert space of square integrable 
complex-valued column "functions") with domain H 1 (M;C m ) (Sobolev space of 
complex- valued column "functions" which are square integrable together with their 
first partial derivatives) and the spectrum of A is discrete. These facts are easily 
established by constructing the parametrix (approximate inverse) of the operator 
A + il. Note that for the special case of the massless Dirac operator a detailed 
examination of relevant functional analytic properties was performed in Chapter 4 

of peg. 

Let Xk and Vk — (vki(x) ■ ■ ■ Vkm( x )) be the eigenvalues and eigenfunctions 
of the operator A. The eigenvalues are enumerated in increasing order with 
account of multiplicity, using a positive index fc = 1,2,... for positive and 
a nonpositive index k = 0,-1,-2,... for nonpositive A^. If the operator A is 
bounded from below (i.e. if mT = 0) then the index fc runs from some integer 
value to +oo; if the operator A is bounded from above (i.e. if rrv = 0) then the 
index k runs from — oo to some integer value; and if the operator A is unbounded 
from above and from below (i.e. if m + ^ and mT ^ 0) then the index k runs 
from — oo to +oo. 
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We will be studying the following three objects. 



Object 1. Our first object of study is the propagator, which is the one- 
parameter family of operators defined as 



U(t):=e- itA = Y,e- itXk v k (x) I [v k (y)]* ( • ) dy , (1.1) 
k Jm 



t e R. The propagator provides a solution to the Cauchy problem 

w \t=o = « 

for the dynamic equation 



(1.2) 



D t w + Aw = 0, (1.3) 

where D t := —id/dt. Namely, it is easy to see that if the column of half-densities 
v = v(x) is infinitely smooth, then, setting w := U(t) v, we get a time-dependent 
column of half-densities w(t, x) which is also infinitely smooth and which satisfies 
the equation (II .3[) and the initial condition (II. 2p . The use of the letter "IT" for the 
propagator is motivated by the fact that for each t the operator U (t) is unitary. 

Object 2. Our second object of study is the spectral function, which is the 
real density defined as 

e(X,x,x):= IM*)H 2 > (1-4) 

0<A fc <A 

where ||i;fc(a;)|| 2 := [ v k(x)]*Vk(x) is the square of the Euclidean norm of the eigen- 
function i>k evaluated at the point x S M and A is a positive parameter (spectral 
parameter). 

Object 3. Our third and final object of study is the counting function 



N(\):= 1 = / e{X,x,x)dx. (1.5) 

0<A fc <A 



In other words, iV(A) is the number of eigenvalues Afe between zero and A. 

It is natural to ask the question: why, in defining the spectral function (|1.4[) and 
the counting function (jl.5l) , did we choose to perform summation over all positive 
eigenvalues up to a given positive A rather than over all negative eigenvalues up to 
a given negative A? There is no particular reason. One case reduces to the other 
by the change of operator A H> —A. This issue will be revisited in Section IT2"1 

Further on we assume that m + > 0, i.e. that the operator A is unbounded 
from above. 



Our objectives are as follows. 
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Objective 1. We aim to construct the propagator (jl.ip explicitly in terms of 
oscillatory integrals, modulo an integral operator with an infinitely smooth, in the 
variables t, x and y, integral kernel. 

Objectives 2 and 3. We aim to derive, under appropriate assumptions on 
Hamiltonian trajectories, two-term asymptotics for the spectral function (|1.4[) and 
the counting function (jl.5[) . i.e. formulae of the type 

e(A, x, x) = a{x) A™ + b(x) A"" 1 + o(A™- 1 ), (1.6) 

N(X) = a\ n + bX' 1 - 1 + o(A"- 1 ), (1.7) 

as A — » +oo. Obviously, here we expect the real constants a, b and real densities 
a(x), b(x) to be related in accordance with 

a = a(x) dx, (1-8) 

J M 

b= I b(x)dx. (1.9) 

It is well known that the above three objectives are closely related: if one 
achieves Objective 1, then Objectives 2 and 3 follow via Fourier Tauberian theo- 
rems [m [571 ED EH- 

We are now in a position to state our main results. 

Result 1. We construct the propagator as a sum of m oscillatory integrals 

U ( t) ™*?~^ u V) ( j t)i (L10 ) 

3 

where the phase function of each oscillatory integral U^(t) is associated with 
the corresponding Hamiltonian h^'(x,^). The symbol of the oscillatory inte- 
gral U^\i) is a complex- valued m x m matrix-function u^\t;y,rf), where y = 
(y 1 , . . . ,y n ) is the position of the source of the wave (i.e. this is the same y that 
appears in formula (ll.ip ) and r\ — . . . , ry„) is the corresponding dual variable 
(covector at the point y). When |?7| — > +oo, the symbol admits an asymptotic 
expansion 

u^(t;y, V ) = 4 j) (t; y, r,) + (t;y, ?,) + ... (1.11) 

into components positively homogeneous in rj, with the subscript indicating degree 
of homogeneity. 

The formula for the principal symbol of the oscillatory integral (t) is known 
|25l |2"U] and reads as follows: 

u^(t;y,r,) = [v^(x^(t;y,r,U^(t;y,r,))] [v^(y, V )}* 

x exp \-iJ q^(x^(r;y,r,),^(T;y,r ] ))dT\ , (1.12) 
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where v^'(z,Q is the normalised eigenvector of the principal symbol A\{z, £) 
corresponding to the eigenvalue (Hamiltonian) hP> (z, £), [x^> (t; y, Tj),^) (f- ^ 77)) 
is the Hamiltonian trajectory originating from the point (y, 77), i.e. solution of the 
system of ordinary differential equations (the dot denotes differentiation in t) 

A® = h { p(x^,^), = -^(i^,^ 1 ) (1-13) 

subject to the initial condition (x (j) , \ t=Q = (y,T)), <? (j) : T'M — > M is the 
function 

q <J) ■- [ v (i)]*A suh v^ - ^{[v U) ]*,A 1 - h u \v u) }-i[v u) ]*{v U \h U) } (1.14) 
and 

A suh (z, C) := A (z, C) + \{A x ) zHa (z, C) (1.15) 

is the subprincipal symbol of the operator A, with the subscripts z a and C, a in- 
dicating partial derivatives and the repeated index a indicating summation over 
a = l,...,n. Curly brackets in formula (|1.14|) denote the Poisson bracket on 
matrix-functions 

:= P z «R Qa - P Ca R z « (1.16) 

and its further generalisation 

{P,Q,R}:=P z «QR Ca -P ia QR z «. (1.17) 

As the derivation of formula (|1.12j) was previously performed only in theses 
|25[ I20|. we repeat it in Sections [2] and [3] of our paper. Our derivation differs 
slightly from that in [25] and [20], 

Formula (|1.12p is invariant under changes of local coordinates on the manifold 
M, i.e. elements of the mxm matrix-function Uq (t; y, rj) are scalars on R x T'M. 
Moreover, formula (|1.12|) is invariant under the transformation of the eigenvector 
of the principal symbol 

V U) ^ e ** W) u O) ) (1.18) 

where 

W : T'M -> K (1.19) 

is an arbitrary smooth function. When some quantity is defined up to the action 
of a certain transformation, theoretical physicists refer to such a transformation as 
a gauge transformation. We follow this tradition. Note that our particular gauge 
transformation (|1.18p . (|1.19p is quite common in quantum mechanics: when <pU) 
is a function of the position variable x only (i.e. when <j>^ : M —> K) this gauge 
transformation is associated with electromagnetism. 

Both Y. Safarov [35] and W.J. Nicoll [5U] assumed that the operator A is semi- 
bounded from below but this assumption is not essential and their formula (|1.12D 
remains true in the more general case that we are dealing with. 

However, knowing the principal symbol (|1.12[) of the oscillatory integral U^\t) 
is not enough if one wants to derive two-term asymptotics fj 1 . 6|) and (|1.7jl . One 
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needs information about u_{(t; y, rf), the component of the symbol of the oscilla- 
tory integral U^\t) which is positively homogeneous in 77 of degree -1, see formula 

(jl.lljl . but here the problem is that y, if) is not a true invariant in the sense 

that it depends on the choice of phase function in the oscillatory integral. We 
overcome this difficulty by observing that U^> (0) is a pseudodifferential operator, 
hence, it has a well-defined subprincipal symbol \U^' (0)] su b- We prove that 

tr[U^(0)] snh = -i{[v^}*,v^} (1.20) 

and subsequently show that information contained in formulae (|1.12p and p. 201) 
is sufficient for the derivation of two-term asymptotics (|1.6p and (|1.7|) . 

Note that the RHS of formula (|1.20l) is invariant under the gauge transforma- 
tion (fTTTgf . ([09]) . 

Formula (|1.20p plays a central role in our paper. Sections [2] and [3] provide 
auxiliary material needed for the proof of formula (|1.20p , whereas the actual proof 
of formula fjl . 20[) is given in Section [4] 

Let us elaborate briefly on the geometric meaning of the RHS of (| 1 .20(1 (a more 
detailed exposition is presented in Section [5J . The eigenvector of the principal 
symbol is defined up to a gauge transformation (|1.18p . (| 1 . 19|) so it is natural to 
introduce a U(l) connection on T'M as follows: when parallel transporting an 
eigenvector of the principal symbol along a curve in T'M we require that the 
derivative of the eigenvector along the curve be orthogonal to the eigenvector 
itself. This is equivalent to the introduction of an (intrinsic) electromagnetic field 
on T'M, with the 2n-component real quantity 

i([«w>]*i£2, [v^Tvg) (i-2i) 

playing the role of the electromagnetic covector potential. Our quantity (jl.21|) is 
a 1-form on T'M, rather than on M itself as is the case in "traditional" electro- 
magnetism. The above U(l) connection generates curvature which is a 2-form on 
T'M, an analogue of the electromagnetic tensor. Out of this curvature 2-form one 
can construct, by contraction of indices, a real scalar. This scalar curvature is the 
expression appearing in the RHS of formula (|1.20|) . 

Observe now that Ylj U^(0) is the identity operator on half-densities. The 
subprincipal symbol of the identity operator is zero, so formula f| 1 . 20[) implies 

2{[«W]*,«W)}=o. (1.22) 

3 

One can check the identity (|1.22l) directly, without constructing the oscillatory 
integrals U^\t): it follows from the fact that the v^(x, £) form an orthonormal 
basis, see end of Section [5] for details. We mentioned the identity (| 1 .22(1 in order 
to highlight, once again, the fact that the curvature effects we have identified are 
specific to systems and do not have an analogue in the scalar case. 

Results 2 and 3. We prove, under appropriate assumptions on Hamiltonian 
trajectories (see Theorems 18.31 and 18.41 for details), asymptotic formulae (|1.6jl and 
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(fLT|) with 



a(x)=Y t J <lt, (1.23) 
j=1 htt)( x ,£)< 



<1 



b(x) = -nJ2 J ([v (j) }*A anh r 
j=1 W) {x,Cj<i 



0) 



n — 1 / 



and a and & expressed via the above densities (|1.23l) and (|1.24[) as (|1.8|l and (|1.9p . 
In (jT^j) and (fLMjl ^ is shorthand for d£ := (27r)-"d£ = (27r)~ n <2£i . . . 
and the Poisson bracket on matrix-functions { • , ■ } and its further generalisation 
{ • , • , • } are defined by formulae (|1.16[) and (11.1 T[) respectively. 

To our knowledge, formula (|1.24p is a new result. Note that in [35] this for- 
mula (more precisely, its integrated over M version (|1.9p ) was written incorrectly, 
without the curvature terms — ;fzy / h^' {[v^']* , See also Section [TBI where 

we give a more detailed bibliographic review. 

It is easy to see that the right-hand sides of (|1.23[) and (jl.24l) behave as densities 
under changes of local coordinates on the manifold M and that these expressions 
are invariant under gauge transformations (|1.18ll , (|1.19p of the eigenvectors of the 
principal symbol. Moreover, the right-hand sides of (|1 .23|) and ()1.24p are unitarily 
invariant, i.e. invariant under transformations of the operator 

A^RAR*, (1.25) 

where 

R: M -> U(m) (1.26) 

is an arbitrary smooth unitary matrix-function. The fact that the RHS of (|1.24p is 
unitarily invariant is non-trivial: the appropriate calculations are presented in Sec- 
tion[Hl The observation that without the curvature terms —^zi f h^{[v^']*, v^'} 
(as in [5S]) the RHS of (|1.24p is not unitarily invariant was a major motivating 
factor in the writing of this paper. 

We will now start making additional assumptions which will, in the end, allow 
us to provide a simple spectral theoretic characterisation of the massless Dirac 
operator. 

Additional assumption 1: 

m = 2 and tr A x = 0. (1.27) 

In this case we can simplify notation by denoting the positive eigenvalue of 

fv + \ 

the principal symbol by h + , the corresponding eigenvector by v + = I \ \ and 



> v. 



2 
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Hamiltonian trajectories by (x + (t; y , 77) , £ + (i; y , 77)) . Obviously, the other eigen- 
value of the principal symbol is — h + , the corresponding eigenvector is 



J i 

and Hamiltonian trajectories are (x + (—t;y,r]),^ + (—t;y,r])) (time reversal). Note 
that in theoretical physics the antilinear transformation 



is referred to as charge conjugation [6]. 

Moreover, in this case the two scalar invariants, — h + ,v + } and 

h + {[v + ]* ,v + }, appearing in formula (|1.24p cease being independent and become 
related as {[«+]*, A x - h+,v+} = ~2h + {[v + ]* ,v+}. Hence, formulae (fl~2"3"f and 
(|1.24p simplify and now read 



a(x) = j <?f, (1.29) 

h+(x,(.)<l 

b(x) = -n J ^[ v +]*A suh v+ + -^jih+{[v+}*,v + }yx,0^. (1.30) 

h+(x,£)<l 

Additional assumption 2: 

the operator A is differential. (1-31) 
In this case there are three further simplifications. 

Firstly the dimension of the manifold can only be n = 2 or n = 3. This follows 
from the ellipticity condition and the fact that the dimension of the real vector 
space of trace- free Hermitian 2x2 matrices is 3. 

Secondly, the subprincipal symbol A sn b does not depend on the dual variable £ 
(momentum) and is a function of x (position) only. 

Thirdly, we acquire a geometric object, the metric. Indeed, the determinant of 
the principal symbol is a negative definite quadratic form 

detA 1 (x,0=-9 aP ^p (1.32) 

and the coefficients g a ^{x), a, /3 = 1, . . . , n, appearing in (|1.32p can be interpreted 
as the components of a (contravariant) Riemannian metric. This implies, in par- 
ticular, that our Hamiltonian (positive eigenvalue of the principal symbol) takes 
the form 

h + (x,€) = y/g°P(x)Z a Z( l (1.33) 

and the ^-components of our Hamiltonian trajectories become geodesies. More- 
over, formulae (jl.29p and (|1.8p simplify and now read 



a(x) = (2Tr)- n u; n ^detg a p(x), (1.34) 
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a = {2n)- n uj n Vo\M , (1.35) 

where u> n is the volume of the unit ball in K" and VolM is the n-dimensional 
volume of the Riemannian manifold M. 

Additional assumption 3: 

u = 3. (1.36) 
In this case there are three more simplifications. 

Firstly, the manifold M is bound to be parallelizable (and, hence, orientable). 
The relevant argument is presented in the beginning of Section [TU] From this point 
we work only in local coordinates with prescribed orientation. 

Secondly, we acquire the identity 

det^ = -iMiAJ^AiMMh,)] 3 (1-37) 
which allows us to define the topological invariant 

c := -~y/det gafl tr^Ai)^)^^,). (1.38) 

The number c defined by formula (|1.38[) can take only two values, +1 or —1, and 
describes the orientation of the principal symbol Ai(x, £) relative to the chosen 
orientation of local coordinates, see formula (|10.4I1 for a more natural geometric 
definition. In calling the number c a topological invariant we are referring to the 
topology of deformations of the elliptic trace- free principal symbol A\(x, £) rather 
than the deformations of the manifold M itself. 

Thirdly, we acquire a new differential geometric object, namely, a teleparallel 
connection. This is an affine connection defined as follows. Suppose we have a 
covector n based at the point y € M and we want to construct a parallel covector 
£ based at the point x € M. This is done by solving the linear system of equations 

A 1 (x,0=A 1 (y,r l ). (1.39) 

Equation (|1.39[) is equivalent to a system of three real linear algebraic equations 
for the three real unknowns, components of the covector £, and it is easy to see 
that this system has a unique solution. It is also easy to see that the affine con- 
nection defined by formula (jl.39l) preserves the Riemannian norm of covectors, 
i.e. g a @(x) ^ a £,/3 — 9 a ^{y) r la r ll3i hence, it is metric compatible. The parallel trans- 
port defined by formula (|1.39p does not depend on the curve along which we trans- 
port the (co)vector, so our connection has zero curvature. The word "teleparallel" 
(parallel at a distance) is used in theoretical physics to describe metric compatible 
affine connections with zero curvature. This terminology goes back to the works 
of A. Einstein and E. Cartan |331 28, 9j, though Cartan preferred to use the term 
"absolute parallelism" rather than "teleparallelism". 

The teleparallel connection coefficients T a p 1 (x) can be written down explicitly 
in terms of the principal symbol, see formula f)10.7[> . and this allows us to define 
yet another geometric object — the torsion tensor 

TV := TV - T a lP . (1.40) 
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Further on we raise and lower indices of the torsion tensor using the metric. Torsion 
is a rank three tensor antisymmetric in the last two indices. Because we are working 
in dimension three, it is convenient, as in [5], to apply the Hodge star in the last 
two indices and deal with the rank two tensor 

T a /3 := X - T^ 5 e l5p y^U^ (1.41) 

rather than with the rank three tensor T. Here e is the totally antisymmetric 
quantity, £123 := + 1. 

The teleparallel connection is a simpler geometric object than the U(l) connec- 
tion because the coefficients of the teleparallel connection do not depend on the 
dual variable (momentum), i.e. they are "functions" on the base manifold M. The 
relationship between the two connections is established in Section [10] where we 
show that the scalar curvature of the U(l) connection is expressed via the torsion 
of the teleparallel connection and the metric as 

- i{[ V +];v+Hx,t) = c - ■ a- 42 ) 

2 [g^ix)^^) 6 ' 2 

Integration of both terms appearing in formula (11.30)) can now be carried out 
explicitly, giving 

([v+]*A suh v+)(x, $t = — L ( tr A sub ^/det g a p ) (x) , (1.43) 

h+(x,£)<l 

h+{[v+]*,v+}(x,Od^^(trTy/d^)(x), (1.44) 
h+(x,£)<l 

where txT := T a a . Note that trT corresponds to one of the three irreducible 
pieces of torsion, namely, the piece which is labelled by theoretical physicists by 
the adjective "axial", see [5J [T3] for details; it is interesting that this is exactly the 
irreducible piece of torsion which is used when one models the neutrino [TU] or 
the electron [8] by means of Cosserat elasticity. Formula (|1.44| follows immedi- 
ately from (| 1 .42(1 . whereas formula (| 1 .43(1 is somewhat less obvious. In order to 
see where formula (|1.43|) comes from one has to write the orthogonal projection 
v+(x,£) [v+(x,Z)}* as v+(x,0 [v + (x,0\* = ^^T) ( A i (*> + h+ 7 ) and use 
the fact that the principal symbol A\(x,£) is an odd function of ^. 
Substituting (fl~36)) . pT43l) and (f04| into (fl~30| we get 

1 * 

= ^2 ( (c tr T - 2 tr A suh ) ^det g af) ){x) . (1.45) 

* 

An explicit self-contained expression for tr T is given in formula (|10.28p . 

Note that the two traces appearing in formula (|1.45j) have a different meaning: 
tr T is the trace of a 3 x 3 tensor, whereas tr A su h is the trace of a 2 x 2 matrix. 
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We now turn our attention to the massless Dirac operator. This operator is 
defined in Appendix see formula (|A.3[) , and it does not fit into our scheme 
because this is an operator acting on a 2-component complex-valued spinor (Weyl 
spinor) rather than a pair of complex-valued half-densities. However, on a par- 
allelizable manifold components of a spinor can be identified with half-densities. 
We call the resulting operator the massless Dirac operator on half- densities. The 
explicit formula for the massless Dirac operator on half-densities is (|A.30[) . 

The massless Dirac operator on half-densities is an operator of the type de- 
scribed in this section (elliptic self-adjoint first order operator acting on a column 
of complex- valued half-densities) which, moreover, satisfies the additional assump- 
tions (|1.27p . (|1.3ip and (| 1 .30[) . We address the question: is a given operator A a 
massless Dirac operator? The answer is given by the following theorem which we 
prove in Section [TT] 

Theorem 1.1 Let A be an elliptic self-adjoint first order pseudodifferential op- 
erator acting on columns of m complex-valued half- densities over a compact n- 
dimensional manifold. Suppose also that this operator satisfies the additional as- 
sumptions 11.27]) . 11.31]) and il.36]) . Then A is a massless Dirac operator on 
half- densities if and only if the following two conditions are satisfied at every point 
of the manifold M : a) the subprincipal symbol of the operator, A su ^(x), is propor- 
tional to the identity matrix and b) the second asymptotic coefficient of the spectral 
function, b{x), is zero. 

The theorem stated above warrants the following remarks. 

• In stating Theorem 11.11 we did not make any assumptions on Hamiltonian 
trajectories (loops). The second asymptotic coefficient (|1.45[) is, in itself, 
well-defined irrespective of how many loops we have. If one wishes to re- 
formulate the asymptotic formula (|1.6|) in such a way that it remains valid 
without assumptions on the number of loops, this can easily be achieved, 
say, by taking a convolution with a function from Schwartz space <S(R). See 
Theorem 1 7. II for details. 

• Conditions a) and b) in Theorem 11.11 are invariant under special unitary 
transformations, i.e. transformations of the operator (| 1 .25[) where R = R(x) 
is an arbitrary smooth special unitary matrix-function. This is not surpris- 
ing as the massless Dirac operator is designed around the concept of SU(2) 
invariance, see Property 4 in Appendix lAl 

• Condition b) in Theorem 11.11 is actually invariant under the action of a 
broader group: the unitary matrix-function appearing in formula Q1.25P does 
not have to be special. 
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2 Algorithm for the construction of the propagator 

We construct the propagator as a sum of m oscillatory integrals (| 1 . 10[) where each 
integral is of the form 

uV)(t) = J e^^^u^H^y^^i^x^^d^^x^^i^dydr,. (2.1) 

Here we use notation from the book [57], only adapted to systems. Namely, the 
expressions appearing in formula (|2.1[) have the following meaning. 

• The function ipV> is a phase function, i.e. a function R x M x T'M — > C 
positively homogeneous in rj of degree 1 and satisfying the conditions 

<p® (t, x; y, rj) = (x- (t; y, rijf (t; y, rj) +0(\x- x& (t; y, rj) | 2 ), (2.2) 

Im tp {j) (t, x; y, rj) > 0, (2.3) 

detcp% fj (t,xW(t;y ) r ) );y,r ) ) J LO. (2.4) 

Recall that according to Corollary 2.4.5 from [27] we are guaranteed to have 
(|2.4p if we choose a phase function 

<p® (t, x; y, V) = (x~ (t; y, rj)) a (i; y, V ) 

+ \c% (*? V, V) (x - (t; y, V )T (x ~ (t; y, rj)f 

+ OQx-x^(t;y, V )\ 3 ) (2.5) 

with complex-valued symmetric matrix-function C^g satisfying the strict 
inequality Im > (our original requirement (|2.3[) implies only the non- 
strict inequality ImC® > 0). Note that even though the matrix- function 
C^p is not a tensor, the inequalities Im > and Im > are invariant 
under transformations of local coordinates x; see Remark 2.4.9 in [27] for 
details. 

• The quantity u^> is the symbol of our oscillatory integral, i.e. a complex- 
valued m x m matrix-function K x T'M — > C TO which admits the asymptotic 
expansion p. lip . The symbol is the unknown quantity in our construction. 

• The quantity d^j) is defined in accordance with formula (2.2.4) from [37] as 

d^&x^rj) := (detV&jV* = \detS m \ 1/2 e i - s(det2 ^ )/i . (2.6) 

Note that in view of (|2.4p our eLyj is well-defined and smooth for x close 
to x^' (t;y,rj). It is known |27| that under coordinate transformations gLo) 
behaves as a half-density in x and as a half-density to the power —1 in y. 

In formula (|2.6p we wrote (det 2 ^^.) 1 / 4 rather than (det y^ 3 2^ ) 1//2 in order 
to make this expression truly invariant under coordinate transformations. 
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Recall that local coordinates x and y are chosen independently and that r\ is 
a covector based at the point y. Consequently, det fx^m changes sign under 

2 (?) 

inversions of local coordinates x or y, whereas det ip)j.i^ retains sign under 
inversions. 

The choice of (smooth) branch of arg(det 2 </3^2 7) ) is assumed to be fixed. 
Thus, for a given phase function tpV> formula (|2.6[) defines the quantity 
duj) uniquely up to a factor e ifc7r / 2 , k — 0,1,2,3. Observe now that if 
we set t — and choose the same local coordinates for x and y, we get 
(p^a (p,y;y,r)) — I. This implies that we can fully specify the choice of 

branch of arg(det 2 <^f2,j ) by requiring that d^u) (0, y\ y, Tj) = 1. 

The purpose of the introduction of the factor d v u) in (|2.ip is twofold. 

(a) It ensures that the symbol is a function on M x T'M in the full 
differential geometric sense of the word, i.e. that it is invariant under 
transformations of local coordinates x and y. 

(b) It ensures that the principal symbol does not depend on the choice 
of phase function <pv) , See Remark 2.2.8 in [27] for more details. 

• The quantity <;V> is a smooth cut-off function IxMx T'M — ► K satisfying 
the following conditions. 

(a) c w) f. »») = on the set {(t,x; \h^(y,T])\ < 1/2}. 

(b) qu) (/j j x; y, rf) = 1 on the intersection of a small conic neighbourhood of 
the set 

{(t,x;y, V ): x = x^(t;y,r))} (2.7) 
with the set {(t,x;y,rj) : \h^(y,rf)\ > 1}. 

(c) sW(t,x',v,Mj) = sW(t,x;y,rj) for \h^[y,rj)\ > 1, A > 1. 

• It is known (see Section 2.3 in [27] for details) that Hamiltonian trajectories 
generated by a Hamiltonian (x, £) positively homogeneous in £ of degree 1 
satisfy the identity 

(a:«r^«=0, (2.8) 

where (4 j) ) Q/3 : = d(xW) a /dr)p. Formulae (J2HD and ([Ml) imply 

l p^(t,x^(t;y,r,y,y,r,)=0. (2.9) 

This allows us to apply the stationary phase method in the neighbourhood 
of the set (|2.7[) and disregard what happens away from it. 



Our task now is to construct the symbols Uq^ (t;y,rj), j — 1, . . . , m, so that our 
oscillatory integrals E/w(t), j = 1, . . . , m, satisfy the dynamic equations 

(A + A(x,D x ))U^(t) mo i C °° (2.10) 
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and initial condition 

J2uW(0) mod =°° I, (2.11) 

3 

where / is the identity operator on half-densities; compare with formulae (jl.3p . 
(|1.2p and Ql.lOp . Note that the pseudodifferential operator A in formula (|2.10p acts 
on the oscillatory integral U(t) in the variable x; say, if A is a differential operator 
this means that in order to evaluate AU^\t) one has to perform the appropriate 
differentiations of the oscillatory integral (|2.ip in the variable x. Following the 
conventions of Section 3.3 of [27], we emphasise the fact that the pseudodifferential 
operator A in formula (|2.10p acts on the oscillatory integral U(t) in the variable x 
by writing this pseudodifferential operator as A(x,D x ), where D x a := —id/dx a . 
We examine first the dynamic equation (|2.10p . We have 

(D t + A(x,D x ))U^(t)=F^(t), 

where F^> (t) is the oscillatory integral 

( t ) = J e*" 5 f& (t, x; y, r?) f « (i, x; y, r?) d y w (t, x;y, V )(-) dy d V 

whose matrix-valued amplitude f^> is given by the formula 

f U) = DtU U) + (^0) + (d (pW )- 1 (D t d (pW ) + u« (2.12) 

where the matrix-function sw (t, x; y, 77) is defined as 

S U) = e -^ U) (d vU) r 1 A(x,D x )(e^ U) d v(j) ). (2.13) 

Theorem 18.1 from |30| gives us the following explicit asymptotic (in inverse 
powers of 77) formula for the matrix-function (|2.13p : 

= (d^r'Y.^^^^i.Dfx^)]^ , (2.14) 

where 

X W (t,z,x;y,T,) = e i ^ t ^^d v u)(t,z ] y,T,), (2.15) 

Tp (j) (t,z,x;y,v) = <P ij) (t,z;y,r)) - (pW (t, x; y, rf) - ^ (i, x; y, 77) (z - x) 13 . (2.16) 
In formula (|2.14jl 

• a := (ai,...,a n ) is a multi-index (note the bold font which we use to 
distinguish multi-indices and individual indices), ot\ := a\\ ■ ■ ■ a n \ , D" := 

D"l ■ ■ ■ D™£, D z p := —id/dz?, 

• A(x, £) is the full symbol of the pseudodifferential operator A written in local 
coordinates x, 

. AW(x,S) := d?A(x,$, 3f := S£ ■ ■ ■ and d u := d/d£p . 
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When \rj\ — > +00 the matrix- valued amplitude f^\t, x; y, 77) defined by formula 
(12. 12|) admits an asymptotic expansion 

(t, x; y, n) = f[ j) (t, x- y, r,) + (t, x; y, 77) + /«> (t, x;y, V ) + ... (2.17) 

into components positively homogeneous in 77, with the subscript indicating degree 
of homogeneity. Note the following differences between formulae (jl.lip and (|2.17p . 

• The leading term in (|2.17[) has degree of homogeneity 1, rather than as in 
(II. lip . In fact, the leading term in (I2.17P can be easily written out explicitly 

f[ j) {t, x; y, 77) = fa® (t, x- y, rj) + A x (x, <pV) ( t) X] y , u « ( t; ^ v ) , (2.18) 

where Ai(x,^) is the (matrix- valued) principal symbol of the pseudodiffer- 
ential operator A. 

• Unlike the symbol u^' (t; y , 77) , the amplitude f^(t,x;y,rj) depends on x. 

We now need to exclude the dependence on x from the amplitude (t, x; y, 77). 
This can be done by means of the algorithm described in subsection 2.7.3 of [27J. 
We outline this algorithm below. 

Working in local coordinates, define the matrix-function ipxr/ in accordance 

with ((pxn)a^ '■= fx"riB an< i then define its inverse (y4"§) _1 from the identity 

{ L P^)a^[{ l Pxri)^ 1 ]i} 1 '■— 5a 1 . Define the "scalar" first order linear differential oper- 
ators 

4 J) V(0/aA a=l,...,n. (2.19) 

Note that the coefficients of these differential operators are functions of the position 
variable x and the dual variable £. It is known, see part 2 of Appendix E in [27], 
that the operators (|2.19p commute: La — L^La, ot,P — 1, . . . ,n. 

Denote L { £ := (L«)^ • • • (L$) a «, {-^ ] ) a := ■ ■ • (-^) a ", 

and, given an r £ N, define the "scalar" linear differential operator 

veL,«w£<*» (i + jr^jjjj^i*) Lf < <2 - 20) 

where \a\ := ai + . . . + a n and the repeated index /3 indicates summation over 
/3 = l,...,n. 

Recall Definition 2.7.8 from [27]: the linear operator L is said to be positively 
homogeneous in 77 of degree p £ K if for any q £ K and any function / positively 
homogeneous in 77 of degree q the function Lf is positively homogeneous in 77 of 
degree p + q. It is easy to see that the operator (I2.20p is positively homogeneous 

(i) 

in 77 of degree —1 and the first subscript in r emphasises this fact. 
Let &^ be the (linear) operator of restriction to x — (t;y, 77), 

'So''' : = i')\x = xU)(t-y,T 1 ) ! ( 2 - 21 ) 
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and let 



e U) r ■.= e ( j) (^Lr (2.22) 



( i) 

for r = 1,2,.... Observe that our linear operators &l r , r — 0, 1, 2, . . ., are pos- 
itively homogeneous in rj of degree — r. This observation allows us to define the 
linear operator 

+ OG 

6« := &U l > (2-23) 

r=0 

where the series is understood as an asymptotic series in inverse powers of r/. 

According to subsection 2.7.3 of [37]; the dynamic equation (|2.10p can now be 
rewritten in the equivalent form 

6W)/0')= o, (2.24) 

where the equality is understood in the asymptotic sense, as an asymptotic expan- 
sion in inverse powers of rj. Recall that the matrix- valued amplitude /00 (t, x\ y, rj) 
appearing in (j!T2"4"| is defined by formulae ([2~T2]) - ([2TT51) . 

Substituting (|2.23[) and f|2 . 1 7|) into (|2.24p we obtain a hierarchy of equations 

e^f^ = o, (2.25) 
6 W/(i) + 6 U)/W =0) (226) 

6^+6^+60^ = 0, 



positively homogeneous in 7/ of degree 1, 0, —1, .... These are the transport equa- 
tions for the determination of the unknown homogeneous components (t; y, 77), 
u_i(t; y, rj), u^ 2 {t\ Vi r ?); ■ • •! °f the symbol of the oscillatory integral (|2.1D . 

Let us now examine the initial condition (|2.11[) . Each operator U^'(Q) is a 
pseudodifferential operator, only written in a slightly nonstandard form. The 
issues here are as follows. 

• We use the invariantly defined phase function (^00 (0, x\ y, rj) = (x — y) a r\ a + 
OQx — y\ 2 ) rather than the linear phase function [x — y) a r\ a written in local 
coordinates. 

• When defining the (full) symbol of the operator U^\t) we excluded the 
variable x from the amplitude rather than the variable y. Note that when 
dealing with pseudodifferential operators it is customary to exclude the vari- 
able y from the amplitude; exclusion of the variable x gives the dual symbol 
of a pseudodifferential operator, see subsection 2.1.3 in [37]. Thus, at t = 0, 
our symbol it 00 (0; y, 77) resembles the dual symbol of a pseudodifferential 
operator rather than the "normal" symbol. 



We have the extra factor d v u) (0, x; y, rj) in our representation of the operator 
t/00(0) as an oscillatory integral. 



17 



The (full) dual symbol of the pseudodifferential operator LTW (0) can be calcu- 
lated in local coordinates in accordance with the following formula which addresses 
the issues highlighted above: 

E { ~^- ( D " d " uU) (o; * v) eiuW {x ' ,v ' n) d ^ (°> x ' y> ♦/)) L=, . ( 2 - 27 ) 

where a;^ (x; y, 77) = <^( J ) (0, x; y, 77) — (x — y)@ rjp . Formula (|2.27p is a version of 
the formula from subsection 2.1.3 of [57], only with the extra factor (— l)\ a K The 
latter is needed because we are writing down the dual symbol of the pseudodiffer- 
ential operator U^> (0) (no dependence on x) rather than its "normal" symbol (no 
dependence on y). 

The initial condition p. lip can now be rewritten in explicit form as 

E E (jD " d " uU) (0; y ' v) e " <3) {x ' y ' v) dvU) (0 ' x; Vi v) ^=y = 1 ' (2 - 28) 

j a 

where / is the m x m identity matrix. Condition (|2.28p can be decomposed into 
components positively homogeneous in rj of degree 0, — 1, —2, . . ., giving us a hierar- 
chy of initial conditions. The leading (of degree of homogeneity 0) initial condition 
reads 

E«o°(o; y, v) = i, (2.29) 

i 

whereas lower order initial conditions are more complicated and depend on the 
choice of our phase functions ipV'. 

3 Leading transport equations 

Formulae (j2~2"Tj) . (|2~T5j) . (l2T2"l) . fTTg]) and the identity £ a h^(x,£) = h^(x,£) 
(consequence of the fact that hV' (x, £) is positively homogeneous in £ of degree 1) 
give us the following explicit representation for the leading transport equation 

[A, ( S W (t; y, 77),C G) (*; V, ri))-h® (x& (t; y, r,),^ (*; », r?))] u« (t; y, r,) = 0. (3.1) 

Here, of course, h^' (x^'(t;y,rf),^'(t;y,rf)) = h^\y,rf). 
Equation Q3.ip implies that 

v$ (t; y, V ) = u« (z« (t; y, 77), ( t; y> 7?) ) [u ,Ci) ( t . yj ^ (3 2) 

where v^'(z, C) is the normalised eigenvector of the principal symbol Ai(z,£) 
corresponding to the eigenvalue h^'(z, £) and iu^ : K x T'M — » C m is a column- 
function, positively homogeneous in r\ of degree 0, that remains to be found. 
Formulae (|2.29p and (|3.2p imply the following initial condition for the unknown 

column- function vr^> : 

w^(0;y,r ] ) = v^(y,i 1 ). (3.3) 
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We now consider the next transport equation in our hierarchy, equation (|2.26p . 
We will write down the two terms appearing in (|2 . 26[) separately. 

In view of formulae (|2~Tgj) and (|2^D| - (|2~2^1) . the first term in (l2~2l)l) reads 



(3.4) 



where we dropped, for the sake of brevity, the arguments (<; y, rj) in Uq and 
x^\ and the arguments {t,x;y,rj) in (p^\ <fx\ frP an d eLo) . Recall that the 
differential operators La are defined in accordance with formula (|2.19l) and the 
coefficients of these operators depend on (t, x; y, rf). 

In view of formulae (|2TT2|) - (|2TT7|1 and (12T2T1) . the second term in (l2~2l))) reads 



(=0) f U) 



(i) 



(d„ )) _1 (A + (-Ai)^^-)^ + A 



,00 



x=xU) 

+ [A 1 -h®]uV\, (3.5) 



where 



(3.6) 



is the matrix-function from (|2.5[) . In formulae (13.51) and (|3.6p we dropped, for the 
sake of brevity, the arguments (t;y,rj) in it^, w'fj, C^2 and x^\ the arguments 
(x^\t;y, r]),^(t;y, ij)) in A , A x , (Ai)^ a , (Ai)^ p and /iW, and the arguments 

(t, x; y, T]) in d^j) and (/^j^ . 

Looking at (|3.4p and (|3 . 5|) we see that the transport equation (|2.26p has a 
complicated structure. Hence, in this section we choose not to perform the analysis 
of the full equation ()2.26|) and analyse only one particular subequation of this 
equation. Namely, observe that equation (|2.26ll is equivalent to m subequations 



[6 uyu) 



6 W /o W ] 



0, 



1*3, 



(3.7) 
(3.8) 



where we dropped, for the sake of brevity, the arguments (x^ (t; y, 77), ft?') (t; y, 77)) 
in lyv)] and [v^] . In the remainder of this section we analyse (sub)equation 
(|3~7f only. 

Equation ()3.7p is simpler than each of the m— 1 equations (|3 . 8[) for the following 
two reasons. 

• Firstly, the term [Ai — h^]u^\ from (13. 5p vanishes after multiplication by 
[oW] * from the left. Hence, equation (|3.7p does not contain u"l\. 
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• Secondly, if we substitute p. 21) into (|3.7[) . then the term with 
vanishes. This follows from the fact that the scalar function 

has a second order zero, in the variable x, at x = x^(t;y,rf). Indeed, we 
have 
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9 ^b (5 T(*'r+^(^i j) ))« w 

_ (#+A 1 (x,^))) ;e 



V U) 



= [«W]'(Ax)..i,W) - fcW + C<J([«W]*(A 1 ) e ,t,W - ftg>) =0, 

where in the last two lines we dropped, for the sake of brevity, the argu- 
ments (a^(i;y,?7),£^(t;y,?7)) in (Ai)^, (Ai)^, /i^2, /i^, and the argu- 
ment (t; y, rf) in cfl (the latter is the matrix-function from formulae (]2.5p 
and (|3.6p ). Throughout the above argument we used the fact that our \v^] * 
and do not depend on x: their argument is (x^'(t;y, rf),^' (t;y, rf)). 



Substituting ([23), (GE3) and (g^J) into we get 

(A+P (j) (t;y,'7))M j) (i;y^)] T = o. 

where 



(3.9) 



p (3) = i [ v 00] 



^ (l-^a W ) (4 j) (#+^(->^' ) )))^ ) 



'm-^M^C^v®. (3.10) 



Note that the ordinary differential operator in the LHS of formula (|3.9p is a 
scalar one, i.e. it does not mix up the different components of the column-function 
w^\t] y, rf). The solution of the ordinary differential equation (|3.9p subject to the 
initial condition (|3.3p is 

(t;y f r)) = vW(y,rj)exi>(-ij^ p (l} (r;y,r?) dr^ . (3.11) 



Comparing formulae (|3.2p , (|3.1ip with formula (|1.12p we see that in order to prove 
the latter we need only to establish the scalar identity 



p {J Ht;y,r ] ) = q^(x^(t;y,r 1 ),^(t;y,r 1 )) : 



(3.12) 
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where qv> is the function (|1.14[) . In view of the definitions of the quantities 
and q(i\ see formulae p.lOp and (|1.14[) , and the definition of the subprincipal 
symbol (|1.15p . proving the identity (|3.12p reduces to proving the identity 



-2 



d 
drjp 



+ [^(a-CA ,£(*))]* ((Ax),^ + (^^Cg) vM(xW,tW). (3.13) 



Note that the expressions in the LHS and RHS of (|3.13p have different structure. 
The LHS of ([3~T3| is the generalised Poisson bracket {[v^]*,Ai - h^,v^}, see 
(fTTT7| . evaluated at z = x^{t;y,rf), ( = ^(t;y,rj), whereas the RHS of (pTTB")) 
involves partial derivatives (in rf) of [x^' (t; y : rf),£^' (f- y ; 77)) (Chain Rule). In 
writing (13 . 1 3|) we also dropped, for the sake of brevity, the arguments (t, x; y, rf) 
in ipt , <fx j iptj , d v u) and the coefficients of the differential operators La an d 
L^j\ the arguments (x^\ in /il^, (^4i)a«^ a and (Ai)^^, and the arguments 
(t;y,r,) in and eg 



Before performing the calculations that will establish the identity (|3 . 1 3|) we 
make several observations that will allow us to simplify these calculations consid- 
erably. 

Firstly, our function pW) (t- y y rf) does not depend on the choice of the phase 
function <pV' (t,x;y,rf). Indeed, if p^(t; y, rf) did depend on the choice of phase 
function, then, in view of formulae p.2p and p. lip the principal symbol of our 
oscillatory integral U^\t) would depend on the choice of phase function, which 
would contradict Theorem 2.7.11 from [37]. Here we use the fact that operators 
(t) with different j cannot compensate each other to give an integral operator 
whose integral kernel is infinitely smooth in t, x and y because all our U^\t) 
oscillate in t in a different way: tp^' (t,x^(t;y,r));y,rf) = —h^(y,rf) and we 
assumed the eigenvalues h^\y,rf) of our principal symbol Ax(y, rf) to be simple. 

Secondly, the arguments (free variables) in (|3.13p are (t]y,rj). We fix an ar- 
bitrary point (t; y, fj) E K x T'M and prove formula (|3.13|) at this point. Put 
(C^)a /3 : = d(£^)a/dvp- According to Lemma 2.3.2 from [27] there exists a local 
coordinate system x such that det(£^) Q ^ ^ 0. This opens the way to the use of 
the linear phase function 



<p& (t, x- y, rf) = {x- x& (t; y, r,)) a & (t; y, r,) 



(3.14) 



which will simplify calculations to a great extent. Moreover, we can choose a local 
coordinate system y such that 



(ti j) )J(t-,y,rj) = Sj 



(3.15) 



which will simplify calculations even further. 
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The calculations we are about to perform will make use of the symmetry 

o#>f°(#v = o^r^v ( 3 - 16 ) 

which is an immediate consequence of formula (|2.8[) . Formula (|3.16[) appears as 
formula (2.3.3) in (27] and the accompanying text explains its geometric meaning. 
Note that at the point (t; y, fj) formula (|3.16ll takes the especially simple form 

{x^{t ] y,f,) = {xfY a {t-,y,^ (3.17) 

Our calculations will also involve the quantity <p^ave Vi v) where x := 
xW(t;y,fj). Formulae (j3~Hj) . (1^5]) . (13~T5)) and ([337]) imply 

(3.18) 



U)\<*Pt 



Further on we denote | := £w) (£; ^). 

With account of all the simplifications listed above, we can rewrite formula 
(|3.13l) . which is the identity that we are proving, as 



{[v^]*,A! -h( j \vM}(x,£) 

d 2 



-{A 1 (x,^)-h^(y, V ) 



(x\?>y 



■ dx a drj a 

— [x — x 

^^(x^-h^v) 

— [x — x 



(x,ri) = (x,fj) 



where ««) = = x^ (t; y , fj) , (AiJ.-c - (^i)*<*£ Q (s, I), fc^. 



(x,J7)=(2,?7) 

(3.19) 

0) _ 



^ = ^(5,0. = * W) and = €«)(t;y,»j). 



Note that the last two terms in the RHS of p,19|) originate from the term with 
d v (j) in p,13[) : we used the fact that d^u) does not depend on x and that 



1 ,7 



(3.20) 



Formula (|3 . 20[) is a special case of formula (3.3.21) from |27j . 

Note also that the term —h^(y, if) appearing (twice) in the RHS of (13. 19)) will 



vanish after being acted upon with the differential operators qJLq^ 
because it does not depend on x. 
We have 



and 



d 2 

dx a dx@ 



teCi)i 



d 2 



dx a drj a 



{x,rf) = {x,fi) 



[^^(ii^^-^-^C^)^ 

+ [«»]*((i 1 ) x «-^2)(^+f;«(^)^), (3.21) 
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teCi)] 



d 2 
dx a dxP 



^ (Ai (x, ) - (z - x& )"if$ (x® , )) « W (x (j) , ^ ) 

' (x,ri) = (x,fj) 

= [^')]*04iW^\ (3.22) 



where (Ix)^ = (^(z^), = /£2(z,f), = vg(x,i) and 5$ = 

vjp(x,i). We also have 

= ^-[««]*(i 1 W^W. (3.23) 
Using formulae (13.231) and (|3.17p we can rewrite formula (|3.2ip as 



^-(^(x,^)) - (z - i^r^^^W)) 



(x, J) ) = (S,77) 



^([^^^IrW^) (3.24) 



Substituting flgJH) and (pT2"2"j) into ([3"T§|) we see that all the terms with {x^ ] ) al3 
cancel out and we get 

{[v^]\A l -h^,v (j) }(x,i) = 

- [C W) ]*((ii)^ Q - h$eJvW - 2[u»]*((i 1 ) x <» - ^ J 2)4f . (3.25) 

Thus, the proof of the identity Q3.13P has been reduced to the proof of the iden- 
tity (pT25j) . 

Observe now that formula (|3.25[1 no longer has Hamiltonian trajectories present 
in it. This means that we can drop all the tildes and rewrite Q3.25P as 

{[vWy,Ai-hV\vW} = 

- - h&) «W) - 2[v (j) ]*(Ai - hW) xa vg , (3.26) 

where the arguments are We no longer need to restrict our consideration to 

the particular point (x, £) = [x, £): if we prove (|3.26p for an arbitrary [x, £) e T'M 
we will prove it for a particular (x, £) € T'M. 

The proof of the identity (|3.26p is straightforward. We note that 

- - h^) x «v£ ~ [v U} ]*(A 1 - hW)e a v$l (3.27) 

and substituting (I3.27P into (I3.26P reduce the latter to the form 

{[vW]*,Ai-hV\vW} = 

[v U} ]*(A 1 - fcCfl) t,C0 _ [«(*)]• (Ai - fr 00 )^?' (3.28) 
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But 

(A! - hW) aa = -[v^ViA, ~ ftW)), (3.29) 

[« W ]*(Aa - ft^) fa = -[^f]*^ - h U) ). (3.30) 
Substituting (ET^ and (HOO} into (EQgf we get 

which agrees with the definition of the generalised Poisson bracket (JTTT7J) . 



4 Proof of formula ( 11.201 ) 



In this section we prove formula (|1.20[) . Our approach is as follows. 
We write down explicitly the transport equations (|3.8p at t — 0, i.e. 

[v d)y [6 0] / p-) + 6 a) / 0-)]|^ = 0j ^ (4il) 

We use the same local coordinates for x and y and we assume all our phase 
functions to be linear, i.e. we assume that for each j we have (|3. 14[) . Using linear 
phase functions is justified for small t because we have )a^(0; V-> v) — and, 
hence, det ^i"^ {t, x; y, rj) 7^ for small t. Writing down equations (|4.1[) for linear 
phase functions is much easier than for general phase functions (|2.2p . 

Using linear phase functions has the additional advantage that the initial con- 
dition (|2.28p simplifies and reads now Y^j u V'(0; y, if) = I. In view of (jl.lip . this 
implies, in particular, that 

5>^(0)=0. (4-2) 

3 

Here and further on in this section we drop, for the sake of brevity, the arguments 

(y, v) in u -\- 

Of course, the formula we are proving, formula (|1.20p . does not depend on our 
choice of phase functions. It is just easier to carry out calculations for linear phase 
functions. 

We will show that (|4.ip is a system of complex linear algebraic equations for 
the unknowns u_\(0). The total number of equations (|4.ip is m 2 — m. However, 
for each j and I the LHS of (|4.1I) is a row of m elements, so (|4.ip is, effectively, a 
system of m(m 2 — m) scalar equations. 

Equation (|4.2p is a single matrix equation, so it is, effectively, a system of m 2 
scalar equations. 

Consequently, the system (|4.ip . (|4.2p is, effectively, a system of m 3 scalar equa- 
tions. This is exactly the number of unknown scalar elements in the m matrices 
u^(0). 

In the remainder of this section we write down explicitly the LHS of (|4.ip and 
solve the linear algebraic system (|4.ip . (14. 2p for the unknowns u^\(0). This will 
allow us to prove formula (|1.20p . 
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Before starting explicit calculations we observe that equations (|4.ip can be 
equivalently rewritten as 

P (i) [ 6 «)/W + 6 «)/«)]| t=o = 0, l*j, (4.3) 

where PW :— [yW (y, rj)} [v® (y, rj)}* is the orthogonal projection onto the eigenspace 
corresponding to the (normalised) eigenvector «W (y, rj) of the principal symbol. 
We will deal with (|4.3p rather than with (14. ip . This is simply a matter of conve- 



4.1 Part 1 of the proof of formula ( I1.20B 

Our task in this subsection is to calculate the LHS of (|4.3[) . In our calculations 
we use the explicit formula (jl.121) for the principal symbol Ug\t;y,r)) which was 
proved in Section [21 

At t = formula (l3~4l reads 



<9 2 



(Ax (x, „) - h& (y,ri)-(x- V yh$ (y, V ))P U) (y, v) 



x=y 



which gives us 



U) 



(4.4) 



In the latter formula we dropped, for the sake of brevity, the arguments (y,r]). 
At t = formula (13. 5p reads 



CO 



o 



t=o 



+ [^ 1 -/ l W] M W(0), (4.5) 

where gw) is the function (| 1 . 14|) and we dropped, for the sake of brevity, the 
arguments (y,rj). Note that in writing down (|4. 5|) we used the fact that 



[(d v u)) 1 d t d v(j 



<P u 'i \(t,x;y,n)=(0,y,y,ri) 



compare with formula (|3.20[) . 

Substituting formulae (|4.4|l and (|4.5p into (|4.3|l we get 



where 



(h® - h^)P®u { H(0) + P {l) B ( o 0) = 0, I £ j, 



(4.6) 



B, 



A Q u) - lh% a + l (A 1 ) ya . na ) PW) - ih<£P$ + i^O^pW. (4.7) 



The subscript in Pq^ indicates the degree of homogeneity in rj. 
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4.2 Part 2 of the proof of formula |L2| 

Our task in this subsection is to solve the linear algebraic system (|4.6p . (|4.2p for 
the unknowns u^\(0). 

It is easy to see that the unique solution to the system (|4.6[) , (|4.2p is 

Summation in (|4.8[) is carried out over all I different from j. 

4.3 Part 3 of the proof of formula ( 1TT201) 

Our task in this subsection is to calculate [£A-"(0)] su b- 
We have 

[^0)W=t£l(0)-~Pft a . (4.9) 

Here the sign in front of | is opposite to that in (| 1 . 1 5|) because the way we write 
U^{0) is using the dual symbol. 

Substituting ((4781) and (|4~71) into (|4~9| we get 

[c/ w (o)] sub =-^« Q +5: 1 



2 y "<> ^ /i0') - «0 



+ P«p + i(ii)v^)-P (i) - 42^« + i(^i)*--P^]) 
v P(')A sub P^)+pO)A sub p(') w (J) v gg N 

^ - WO 2 V + 2^ h (j) - wo J ' ^- iUj 



where 



G j7 := pCO^^pW _ 2/iWpW + 2(A 1 ) ( -P 1 W] 

+ p 0) [(A W ^ (0 - 2<o p(0 + 2(^0^ P« ] . 



2G 



We have 



Gji = 2P (l) {A 1 , P w } + 2P (j} {A 1 ,P (l) } 

+ P®[(Ai - h^) yar]a pW + 2{Ai - h^) na P$] 
+ P«[(Ax - h®) yana P® + 2(A X - hW)r, a P$] 
= 2P® {Ai, P^ } + 2P^{Ai,P^} — P^{Ai — h^\P^} — P^{Ai — h^,P^} 
+ pOp! - ^)) y «, a P« + (A, - hW) Va P$ + (A, - hW) ya P$] 
+ P«[(A a - hW) w «, a PW + (Ax - hW) Va P$ + (A 1 - h®) y c.P®] 
= P^{Ai + h li) ,P ii) } + P <J) {A 1 + h w ,P ( - l) } 

- pW(Ax - h^)P$ Va - pW(Ai - hW)PX a 
= P^{Ai + h^\P^} + P<j){Ai + /iW,p(0} 

- P«(/*W - fcW)P$,„ - P«(/iW - /i^P^ 

= pW{A 1+ M J ^p«}+pW{A+/iW,p(0} + (/ l (i)_/ l W)(p(Op«^_p 

so formula (|4.10p can be rewritten as 

P« (2A sub pW + i{ Ai + fcCi) , pO) }) + pO') (2A sub P« + i{ A x + fc« , P« }) 



IE 



2 - fr« 

(4.11) 

But 

J2(pwp y % a -p^p$ v j - (E pW ) p ^ - p0) (E p(0 ) . 

= (I- P®)P$ na - P«(7 _ P«)^ = 
so formula (|4. can be simplified to read 

[£/ (j) (0)]sub 

_ l^P«(2A sub P« +i{A 1 + /i«,P«})+P«(2A sub pW +i{A 1 + h«\P®}) 



2 ^ - fc(Q 



4.4 Part 4 of the proof of formula (11.201) 

Our task in this subsection is to calculate ti[U^ (0)] sub . 



(4.12) 
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Formula (|4.12[) implies 

tr ^ Wl sub = 2 tr ^ ' ( 413 ) 

Put A 1 = J2 k /i (fe) P (fc) and observe that 

• terms with the derivatives of h vanish and 

• the only k which may give nonzero contributions are k — j and k — I. 
Thus, formula (|4.13[) becomes 

x(hW[p®{pW ,p(J)} + pU) {P U) \p«)}] +h W[p®{p®,pW}+pW{pV),p<r)}]). 

(4.14) 

We claim that 

tr (p(0 { P U) f P U) }) = tr (p(i) |p(j) j p(0}) 

= -tr(pW{pW,pW}) = -tr(pW{pW,pW}) 

- ([« ( °]*«?- ) )([«^]*' ,(0 ) ~ ([^l^X^]*^ )- ( 4 - 15 ) 

These facts are established by writing the orthogonal projections in terms of the 
eigenvectors and using, if required, the identities 

[4'2]*u (i) + = o, [«W]*» W) + = °> 

In view of the identities (|4.15|) formula (|4.14D can be rewritten as 
tr[tfW(0)Ub = itr £ pW{pW ; pa) } 

= t tr({pW , P^ } - pCi) {P« , Ptf ) }) = -i tr(P« {P W , P W }). (4.16) 
It remains only to simplify the expression in the RHS of (|4.16[) . We have 
tr (P«{P«,P«}) = 

+ [([w^r^xi^r^) - ([« w) ]*«^ ) )([« w) ]*«?- ) )] 

+ [{[v {j) }*v$)([vW]* v W) - (lv U) ]*v^)(lv$]*v^)] 
= {[««)]•, W C7)} + [([«W]*»yj)([t;W)]*t;W) - ([w^l^Xfv^]*^)] 

= {[«W)]*, W W)}. (4.17) 
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Formulae (|3~T5|l and P~TT]) imply formula (|O0]) , 



5 U(l) connection 

In the preceding Sections [2H4] we presented technical details of the construction of 
the propagator. We saw that the eigenvectors of the principal symbol, v^(x, £), 
play a major role in this construction. As pointed out in Section [TJ each of these 
eigenvectors is defined up to a U(l) gauge transformation (|1.18p . (|1.19p . In the 
end, the full symbols of our oscillatory integrals U^\t) do not depend on 

the choice of gauge for the eigenvectors (x, £). However, the effect of the gauge 
transformations (|1.18[) . (| 1 . 19[) is not as trivial as it may appear at first sight. We 
will show in this section that the gauge transformations (|1.18p , (|1.19p show up, 
in the form of invariantly defined curvature, in the lower order terms u_\(t; y : rf) 
of the symbols of our oscillatory integrals U^'(t). More precisely, we will show 
that the RHS of formula (|1.20[) is the scalar curvature of a connection associated 
with the gauge transformation (|1.18p . (|1.19p . Further on in this section, until the 
very last paragraph, the index j enumerating eigenvalues and eigenvectors of the 
principal symbol is assumed to be fixed. 

Consider a smooth curve T C T'M connecting points (y, rf) and (x, £). We write 
this curve in parametric form as (z(t),((t)), t € [0, 1], so that (z(0),£(0)) = {y,rj) 
and (*(1),C(1)) = (z,0- Put 

w (t) :=e^ ( *V 3) (z(i),C(i)), (5.1) 

where cj> : [0, 1] — > R is an unknown function which is to be determined from the 
condition 

iw*w = (5.2) 

with the dot indicating the derivative with respect to the parameter t. Substituting 
(|5.ip into (|5.2p we get an ordinary differential equation for <j> which is easily solved, 
giving 

0(1) = 0(0) + f (z a (t) P a (z(t),C(t)) + Cr(t) Qpf{z(t), C(t))) dt 
Jo 

= 0(0) + J (P a dz a + Q^C 7 ) . (5-3) 

where 

P a :=i[v®]*v$!, CP :=i[vW]*v$. (5.4) 

Note that the 2n-component real quantity (P a , Q 1 ) is a covector field (1-form) on 
T'M. This quantity already appeared in Section Q] as formula (|1.2ip . 
Put f(y,r)) := e^(°), f(x,£) := e^M and rewrite formula (JOJ) as 

fM - fiy^e^riP^+Q-"^). ( 5. 5) 

Let us identify the group U(l) with the unit circle in the complex plane, i.e. with 
/ G C, ll/H = 1. We see that formulae (15 ,5|) and (|5.4p give us a rule for the parallel 
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transport of elements of the group U(l) along curves in T'M . This is the natural 
U(l) connection generated by the normalised field of columns of complex- valued 
scalars 

v^(z,0 = (v{ j) (z,Q ■■■ v$(z,Q) T . (5.6) 

Recall that the T appearing in formula (|5 . 5[) is a curve connecting points (y, rj) 
and (x, £), whereas the v^\z,C) appearing in formulae (|5.4p and (|5.6[) enters our 
construction as an eigenvector of the principal symbol of our mxm matrix pseudo- 
differential operator A. 

In practice, dealing with a connection is not as convenient as dealing with the 
covariant derivative V. The covariant derivative corresponding to the connection 
(|5.5p is determined as follows. Let us view the (x, £) appearing in formula (|5.5p as 
a variable which takes values close to (y, 77), and suppose that the curve T is a short 
straight (in local coordinates) line segment connecting the point (y, rj) with the 
point (a;, £). We want the covariant derivative of our function f(x, £), evaluated at 
(y,T)), to be zero. Examination of formula (|5.5p shows that the unique covariant 
derivative satisfying this condition is 

V Q :=d/dx a -iP a (x,£), V 7 :=d/d^-iQi{x^). (5.7) 

We define the curvature of our U(l) connection as 

d _ _„■ (VctVp — V/jVq V q V 5 -V 5 V q \ . , 

It may seem that the entries of the (2n) x (2n) matrix (|5.8|) are differential opera- 
tors. They are, in fact, operators of multiplication by "scalar functions". Namely, 
the more explicit form of (|5.8jl is 



n> _ a^P dx a d£s dx<* I (r 
\ dxP d£~, d£s d£~, J 

The (2n) x (2n) - component real quantity (|5.9|) is a rank 2 covariant antisymmetric 
tensor (2-form) on T'M. It is an analogue of the electromagnetic tensor. 

Substituting (|5.4[) into (|5.9p we get an expression for curvature in terms of the 
eigenvector of the principal symbol 

Examination of formula (|5 . 10|) shows that, as expected, curvature is invariant 
under the gauge transformation (|1.18p . (|1.19p . 

It is natural to take the trace of the upper right block in (|5.8p which, in the 
notation (|1.16p . gives us 

-i(V Q V Q -V a V Q ) = -i{[v^]*,v^}. (5.11) 

Thus, we have shown that the RHS of formula f| 1 . 20[) is the scalar curvature of our 
U(l) connection. 
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We end this section by proving, as promised in Section[TJ formula (| 1 . 22[) without 
referring to microfocal analysis. In the following arguments we use our standard 
notation for the orthogonal projections onto the eigenspaces of the principal sym- 
bol, i.e. we write pW := v {k ^[v^}*. We have tr{pW,?0')} = o and £,PW = I 
which implies 

hi 

= J2 tv ( pU) { p(j \ p(j) })+ tr(P (i) {P (j) ,P (j) }). (5.12) 

But, according to formula (|4.15p . for ! ^ j we have 

tr(P«{P«,P«}) = -tr(P«{pW,P«}), 

so formula (|5.12p can be rewritten as ^ ■ tr(Pw)^pC?) ) PW}) — 0. It remains only 
to note that, according to formula (0TT7]), tr(P^{pW, P^}) = {[v^]*,v^}. 

6 Singularity of the propagator at t = 

Following the notation of [27], we denote by 

J^t [/(A)] =/(*)= | e-« A /(A)dA 
the one-dimensional Fourier transform and by 

[/(*)] = fW = i^r 1 J e UX f(t)dt 

its inverse. 

Suppose that we have a Hamiltonian trajectory (x^'(t;y,rf),^'(t;y,rf)) and 
a real number T > such that x^ (T; y, rf) — y. We will say in this case that we 
have a loop of length T originating from the point y £ M. 

Remark 6.1 There is no need to consider loops of negative length T because, 
given a T > 0, we have x^(T; y, n + ) = y for some n + £ T'M if and only if we 
have x^(—T;y,r]~) — y for some rf~ £ TyM. Indeed, it suffices to relate the rf 1 
in accordance with rf^ = £W (±T; y, rj^). 

Denote by c R the set of lengths T > of all possible loops generated 
by the Hamiltonian h^\ Here "all possible" refers to all possible starting points 
(y, rj) £ T'M of Hamiltonian trajectories. It is easy to see that T~V) . We put 




infT w if T (j ' ) ^0, 
+oo if = 0. 
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In the Riemannian case (i.e. the case when the Hamiltonian is a square root of 
a quadratic polynomial in £) it is known [23l [21] that there is a loop originating 
from every point of the manifold M and, moreover, there is an explicit estimate 
from above for the number . We are not aware of similar results for general 
Hamiltonians. 

We also define T := min T^'). 

j — l,...,m+ 

Remark 6.2 Note that negative eigenvalues of the principal symbol, i.e. Hamil- 
tonians hU'(x,£) with negative index j = — 1, . . . , —mT, do not affect the asymp- 
totic formulae we are about to derive. This is because we are dealing with the case 
A — > +oo rather than A — > — oo. 

Denote by 

u(t,x,y) :=J2 e ~ itXkv k( x )l v k(y)}* (6-1) 
k 

the integral kernel of the propagator . The quantity (|6.ip can be understood 
as a distribution in the variable t £ R depending on the parameters x, y G M. 
The main result of this section is the following 

Lemma 6.1 Let p : K — > C be an infinitely smooth function such that 

suppp C (— T,T), (6.2) 

p(0) = 1, (6.3) 
p'(0) = 0. (6.4) 
Then, uniformly over y £ M , we have 

T-_l x [p(t) tr u(t, y, y)} = n a(y) A" -1 + (n - 1) b(y) A"~ 2 + 0(A"~ 3 ) (6.5) 



as A — > +oo. The densities a(y) and b(y) appearing in the RHS of formula i6.5]) 
are defined in accordance with formulae $1.23\) and \1.2J$ . 

Proof Denote by (S*M)^ the (n— l)-dimensional unit cosphere in the cotan- 
gent fibre defined by the equation h^'(y,rf) = 1 and denote by d(S*M)^' the 
surface area element on (S*M)^' defined by the condition dr\ = d(S*M)^ dh^\ 
The latter means that we use spherical coordinates in the cotangent fibre with the 
Hamiltonian ftW) playing the role of the radial coordinate, see subsection 1.1.10 of 
[27] for details. In particular, as explained in subsection 1.1.10 of [27], our surface 
area element d{S*M)^ is expressed via the Euclidean surface area element as 

/ ™ \-i/2 
d(S*M)^ = ^(^(l/,*?)) x Euclidean surface area element . 

Denote also d(S*M)^ := (2%)- n d(5*M)W . 
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According to Corollary 4.1.5 from [37] we have uniformly over y € M 
•^A[p(*) tru (*'2/>y)] = 

m+ 

(c 0) (y) A"" 1 + (y) A"~ 2 + e« (y) A"~ 2 ) + 0(A"- 3 ) , (6.6) 



where 

C W(y)= J trui j \0;y,r))a(S;M)W, (6.7) 

(n-1) J tr(^-z4 j) (0;y,r ? ) + ^{4 j) | t=0 ,/ l «}(y,^^M)W, (6.8) 
e i3 \y)= J tr[U^(0)} sub (y,r,)d(S* y M)^ . (6.9) 

(S*Af)U) 

Here uffi (t; y, rf) is the principal symbol of the oscillatory integral (|2 . 1 1) and ii^ (t; y, rf) 
is its time derivative. Note that in writing the term with the Poisson bracket in 
(|6.8p we took account of the fact that Poisson brackets in [127] and in the current 
paper have opposite signs. 

Observe that the integrands in formulae (|6.7[) and (j6.8[) are positively homo- 
geneous in -q of degree 0, whereas the integrand in formula (j6.9[) is positively 
homogeneous in r\ of degree —1. In order to have the same degree of homogeneity, 
we rewrite formula (|6.9p in equivalent form 



e (j) (2/)= / (h (j hY[U^(0)} suh )(y,f 1 )d(S;M)^ . (6.10) 

Switching from surface integrals to volume integrals with the help of formula 
(1.1.15) from [57], we rewrite formulae IfBTT)) . ([6T5]) and ([BTTU]) as 

c u) (y)=n J trutf ) (Q;y,r))dv, (6.11) 

fc(i) (y,r))<l 

d (j \y) = n(n-l)x 

J tr(-iu^\0;y,r,) + ^{u ( j) \ t=o ,h^}(y,v)^V, (6.12) 

/iU)(y,77)<l 
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e^(y) = n I (h® tr[L/«(0)] sub ) (y, rj) dr) . (6.13) 

h(J) (y,ri)<l 

Substituting formulae (j!.12|l and (|1.14[) into formulae (|6.1ip and (|6.12ll we get 




(6.14) 



d^(y) = - n (n-l)x 

[v^]*A suh v^ - '-{[v^WA^h^M^ (y,r])<tr)- (6-15) 

Substituting formula (|1.20|) into formula (|6.13| we get 

e {l) (y) = -ni J (h^{[v u) \* ,v u} })(y, rf) tir) . (6.16) 

W>(y,r))<l 

Substituting formulae (|6.14[) - (|6.16p into formula (|6.6p we arrive at (|6.5[) . □ 

Remark 6.3 The proof of Lemma \6.1\ given above was based on the use of Corol- 
lary 4-1-5 from ]27$ . In the actual statement of Corollary 4. 1.5 in f27f uniformity 
in y £ M was not mentioned because the authors were dealing with a manifold 
with a boundary. Uniformity reappeared in the subsequent Theorem 4-2.1 which 
involved pseudodifferential cut-offs separating the point y from the boundary. 

7 Mollified spectral asymptotics 

Theorem 7.1 Let p : M — > C be a function from Schwartz space S(M) whose 
Fourier transform p satisfies conditions Ii6.2]) - [6^4\ l. Then, uniformly over x G M , 
we have 

f e (\-V,x,x)p(p)dv = a(x)\ n +b(x)\ n - 1 + r^ n ~ 2) l [ n ~ 3 > (7.1) 
J ^O(lnA) if n = 2, 

as A — > +00. The densities a(x) and b(x) appearing in the RHS of formula i7.1]) 
are defined in accordance with formulae il.23\) and \1.2J$ . 

Proof Our spectral function e(A, x, x) was initially defined only for A > 0, see 
formula (|1.4j) . We extend the definition to the whole real line by setting 

e(A, x, x) := for A < 0. 
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Denote by e'(A, x, x) the derivative, with respect to the spectral parameter, of 
the spectral function. Here "derivative" is understood in the sense of distributions. 
The explicit formula for e'(A, x, x) is 



+00 



e'(A,x,x) :=X> fc (s)|| a *(A-A fc ). (7.2) 



fc=l 

Formula (|7.2[) gives us 

+00 



e'(A - p{p)dp = ^ \\vk(x)\\ 2 p{\ - Afe). (7.3) 

fe=i 

Formula (|7.3p implies, in particular, that, uniformly over x € M , we have 

e'{X- p,x,x)p(p)dp = 0(|A|"°°) as A -» -00 , (7.4) 



where 0(|A| 00 ) is shorthand for "tends to zero faster than any given inverse power 
of|A|". 

Formula (|7.3|) can also be rewritten as 

e'(A - /i, x, x) p(fi) dp = J^ x [p(t) tr u(t, x, x)] - \\v k (x)\\ 2 p(X - X k ) , (7.5) 

fc<0 

where the distribution u(t,x,y) is defined in accordance with formula (|6.1I) . Clearly, 
we have 

^||u fe (x)|| 2 p(A-A fe )=0(A- 00 ) as A -> +00 . (7.6) 

fe<0 

Formulae (|7.5|) . (I7.6[) and Lemma IUTT1 imply that, uniformly over a; € M, we have 

e (A — /i, a;, a;) p(p) dpi = 

na(x)X n - 1 + {n-l)b{x)X n - 2 +0(X n - 3 ) as A — > +00 . (7.7) 
It remains to note that 

— [ e(X- p,x,x) p(p)dp = ( e'{X- p,x,x) p(p)dp. (7.8) 



dX . 

Formulae dTTHJ), {HJ) and ([721) imply ([TIT]) . □ 

Theorem 7.2 Lei p : 1 -> C 6e a function from Schwartz space <S(R) whose 
Fourier transform p satisfies conditions i6.2]) - [K^\ ). Then we have 

N(X~p)p(p)dp = aX n + bX n - 1 + \° {Xn ^ ) lf n ~ 3, (7.9) 



O(lnA) if n = 2, 

as X — > +00 . The constants a and b appearing in the RHS of formula \7.9\) are 
defined in accordance with formulae hi. 8)) , hi. 23]) , lil.9\) and 11.24]) - 
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Proof Formula (|7.9p follows from formula (|7.ip by integration over M, see also 
formula ([L5]) . □ 

In stating Theorems 17.11 and 17.21 we assumed the molliher p to be complex- 
valued. This was done for the sake of generality but may seem unnatural when 
mollifying real- valued functions e(A, x, x) and N(X). One can make our construc- 
tion look more natural by dealing only with real-valued mollifiers p. Note that 
if the function p is real-valued and even then its Fourier transform p is also real- 
valued and even and, moreover, condition (|6.4[) is automatically satisfied. 

8 Unmollified spectral asymptotics 

In this section we derive asymptotic formulae for the spectral function e(X,x,x) 
and the counting function N(X) without mollification. The section is split into 
two subsections: in the first we derive one-term asymptotic formulae and in the 
second — two-term asymptotic formulae. 

8.1 One-term spectral asymptotics 
Theorem 8.1 We have, uniformly over x G M , 

e(\,x,x) = a(x)\ n + 0(X n - 1 ) (8.1) 

as X — > +oo. 

Proof The result in question is an immediate consequence of formulae (|7.8p , 
(|7.7p and Theorem 1 7 . 1 1 from the current paper and Corollary B.2.2 from [27] • □ 

Theorem 8.2 We have 

N(X) = a\ n + 0{X n ' 1 ) (8.2) 

asA-> +oo. 

Proof Formula (|8.2p follows from formula (|8.ip by integration over M, see also 
formula HH). □ 

8.2 Two-term spectral asymptotics 

Up till now, in Section [7] and subsection 18. 1[ our logic was to derive asymptotic 
formulae for the spectral function e(A, x, x) first and then obtain corresponding 
asymptotic formulae for the counting function -/V(A) by integration over M . Such 
an approach will not work for two-term asymptotics because the geometric condi- 
tions required for the existence of two-term asymptotics of e(A, a;, x) and N(X) will 
be different: for e{X,x,x) the appropriate geometric conditions will be formulated 
in terms of loops, whereas for N(X) the appropriate geometric conditions will be 
formulated in terms of periodic trajectories. 

Hence, in this subsection we deal with the spectral function e(X,x,x) and the 
counting function N(X) separately. 
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In what follows the point y £ M is assumed to be fixed. 

Denote by ILy the set of normalised (h^\y, rf) = 1) covectors 77 which serve as 
starting points for loops generated by the Hamiltonian hy>. Here "starting point" 
refers to the starting point of a Hamiltonian trajectory (xV> (t; y , rf) , £^ (t; y , 77)) 
moving forward in time (i > 0), see also Remark 16. II 

The reason we are not interested in large negative t is that the refined Fourier 
Tauberian theorem we will be applying, Theorem B.5.1 from [57], does not require 
information regarding large negative t. And the underlying reason for the latter 
is the fact that the function we are studying, e(A, x, x) (and, later, N(X)), is 
real-valued. The real-valuedness of the function e(\,x,x) implies that its Fourier 
transform, e(t, x, x), possesses the symmetry e(— t,x, x) = e{t,x,x). 

The set II^ is a subset of the (n — l)-dimensional unit cosphere (S*M)^ and 
the latter is equipped with a natural Lebesgue measure, see proof of Lemma |6. 11 
It is known, see Lemma 1.8.2 in [27j . that the set is measurable. 

Definition 8.1 A point y € M is said to be nonfocal if for each j = 1, . . . ,m + 
the set Hy^ has measure zero. 

With regards to the range of the index j in Definition 18.11 as well as in sub- 
sequent Definitions I8.2H8.41 see Remark 16.21 

We call a loop of length T > absolutely focused if the function 

\x^(T;y 1 i 1 )-y\ 2 

has an infinite order zero in the variable 77, and we denote by (Hy)^ the set of 
normalised (h^(y, rf) = 1) covectors 77 which serve as starting points for absolutely 
focused loops generated by the Hamiltonian h^\ It is known, see Lemma 1.8.3 in 
[27], that the set (U^l is measurable and, moreover, the set U y 3) \ (11^)® has 
measure zero. This allows us to reformulate Definition 18.11 as follows. 

Definition 8.2 A point y G M is said to be nonfocal if for each j = 1, . . . ,to + 
the set (II")' 3 ' has measure zero. 

In practical applications it is easier to work with Definition 18.21 because the set 
(n°)( J ) is usually much thinner than the set Tly . 

In order to derive a two-term asymptotic formula for the spectral function 
e(A, x, x) we need the following lemma (compare with Lemma |6.1[) . 

Lemma 8.1 Suppose that the point y € M is nonfocal. Then for any complex- 
valued function 7 € C^°(R) with supply C (0, +00) we have 

^ x mt)tru(t,y,y)} = o(A"- 1 ) (8.3) 

as A — > +00 . 

Proof The result in question is a special case of Theorem 4.4.9 from [57] . □ 
The following theorem is our main result regarding the spectral function e(A, x, x). 
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Theorem 8.3 If the point x £ M is nonfocal then the spectral function e(A, x, x) 
admits the two-term asymptotic expansion (1.6\) as X — > +00. 

Proof The result in question is an immediate consequence of formulae (17. 7p , 
Theorem 17.11 and Lemma 18.11 from the current paper and Theorem B.5.1 from 
[27]. □ 

We now deal with the counting function N(X). 

Suppose that we have a Hamiltonian trajectory (x^(t;y ) rf),^'(t' } y ) rf)) and 
a real number T > such that (x^> (T; y, rj),^' (T; y, rf)) — (y, vf). We will say in 
this case that we have a T-periodic trajectory originating from the point (y, rf) € 
T'M. 

Denote by {S*M)V> the unit cosphere bundle, i.e. the (2n — l)-dimensional 
surface in the cotangent bundle defined by the equation h^\y,rf) = 1. The 
unit cosphere bundle is equipped with a natural Lebesgue measure: the (2n — 1)- 
dimensional surface area element on (S*M)^ is dy d(S*M)^ where d{S*M)^ 
is the (n — l)-dimensional surface area element on the unit cosphere (S*M)^', see 
proof of Lemma 16.11 

Denote by 11^ the set of points in (S*M)^ which serve as starting points for 
periodic trajectories generated by the Hamiltonian h^K It is known, see Lemma 
1.3.4 in [27], that the set II^) is measurable. 

Definition 8.3 We say that the nonperiodicity condition is fulfilled if for each 
j = 1, . . . , m + the set IJy) has measure zero. 

We call a T-periodic trajectory absolutely periodic if the function 

\xU)(T;y,T})-y\ 3 + \zW(T ]y ,Ti)-r,\ 3 

has an infinite order zero in the variables (j/,77), and we denote by (II a )W the 
set of points in (S*M)V' which serve as starting points for absolutely periodic 
trajectories generated by the Hamiltonian h^K It is known, see Corollary 1.3.6 in 
[27], that the set (H a )^ is measurable and, moreover, the set \ (n a )(j) has 
measure zero. This allows us to reformulate Definition 18.31 as follows. 

Definition 8.4 We say that the nonperiodicity condition is fulfilled if for each 
j = 1, . . . ,m + the set (JJ a )^ has measure zero. 

In practical applications it is easier to work with Definition 18.41 because the set 
(n°)w is usually much thinner than the set n^). 

In order to derive a two-term asymptotic formula for the counting function 
N(X) we need the following lemma. 

Lemma 8.2 Suppose that the nonperiodicity condition is fulfilled. Then for any 
complex-valued function 7 6 C^°(K) with supply C (0, +00) we have 

I J r rJ> x m)teu(t,y,y)]dy = o(X n - 1 ) (8.4) 
Jm 

as X — > +00. 
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Proof The result in question is a special case of Theorem 4.4.1 from [57]. □ 
The following theorem is our main result regarding the counting function N(X). 

Theorem 8.4 If the nonperiodicity condition is fulfilled then the counting func- 
tion N(\) admits the two-term asymptotic expansion il.7\) as A — > +oo. 

Proof The result in question is an immediate consequence of formulae (jl.5p , 
(|7.7p . Theorem 17.11 and Lemma [8.21 from the current paper and Theorem B.5.1 
from [27J. □ 

9 U(m) invariance of the second asymptotic coeffi- 
cient 

We prove in this section that the RHS of formula (|1.24l) is invariant under unitary 
transformations (|1.25p , (|1.26p of our operator A. The arguments presented in this 
section bear some similarity to those from Section[SJ the main difference being that 
the unitary matrix-function in question is now a function on the base manifold M 
rather than on T'M . 

Fix a point x € M and an index j (index enumerating the eigenvalues and 
eigenvectors of the principal symbol) and consider the expression 



- -\[v^]*,A 1 + — !_ftCj){[«W)]*,««)}Va;,f)de, (9-1) 

2 n — 1 / 

compare with (|1.24p . We will show that this expression is invariant under the 
transformation ([PS]) , fPS]) , 

The transformation (| 1 . 25|) , p.26p induces the following transformation of the 
principal and subprincipal symbols of the operator A: 

Ax h-> RAiR*, (9.2) 

A suh ^ RA snh R* + % - (R x a(Ai)^ a R* - R(Ax) 6a R* xa ) . (9.3) 

The eigenvalues of the principal symbol remain unchanged, whereas the eigen- 
vectors transform as 

V U)^,R V <J). (9.4) 
Substituting formulae (|9.2 p - (|9.4p into the RHS of (|9.ip we conclude that the in- 
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crement of the expression (|9.1I) is 

J ^ [v U) V{R * Rxa{Alha _ ( Al )z a R* xa R)vW 

\ ([vW]*RZ*R(Ai - h&)vg - - h^)R*R x .v^ 

-h {j) (lv U) }*R:«Rv£ - [v^]*R*R xa v^\x,Odt, 



n — 

which can be rewritten as 



I 

2 

h.O)(x,g)<l 



fcgj ([u^-R^ito^ - [v^}*R*R x a,v^ 



In view of the identity R*R = I the above expression can be further simplified, so 
that it reads now 



J U^[ v W]*R*R x <,vW 

L_ ft W) jj*^-^ + [ v ^]*R*R xaV 0)y\( x ,0 d£ . (9.5) 



h(i~>(x,£)<l 



n 

Denote 



B a (x) := -iR*R x « (9.6) 



and observe that this set of matrices, enumerated by the tensor index a running 
through the values is Hermitian. Denote also b a (x,£) := [v^']*B a v^ 

and observe that these b a are positively homogeneous in £ of degree 0. Then the 
expression (|9.5I) can be rewritten as 



/iW>(:e,£)<1 



(»S».-^T» W ^)M«- <* 



7) 



Lemma 4.1.4 and formula (1.1.15) from [57] tell us that the expression (|9.7|) is 
zero. 



10 Teleparallel connection 

In this section we work under the additional assumptions (jl.27p , (|1.3ip and (|1.36p , 
i.e. we study a 2 x 2 matrix differential operator in dimension 3 with trace-free 
principal symbol. Our aim is to show that in this case the principal symbol 
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generates additional geometric structures which allow us to reformulate the results 
of our spectral analysis in a much clearer geometric language. 

Let us show first that the manifold M in this case is parallelizable. The prin- 
cipal symbol A\ (x, £) is linear in f so it can be written as 

A t {x,0 =<r a (x)£ a , (10.1) 

where a a (x), a — 1,2,3, are some trace-free Hermitian 2x2 matrix-functions. 
Let us denote the elements of the matrices a a as a a ab where the dotted index, 
running through the values 1,2, enumerates the rows and the undotted index, 
running through the values 1,2, enumerates the columns; this notation is taken 
from [TU]. Put 

V 1 a (x):^Rea a i2 (x), V 2 a (x) := - Ima a i2 (x), V 3 a (x) := a a h (x). (10.2) 

Formula (|10.2[) defines a triple of smooth real vector fields Vj(x), j = 1,2,3, on 
the manifold M. These vector fields are linearly independent at every point x 
of the manifold: this follows from the fact that det Ai (x,£) ^ 0, V(x,£) £ T'M 
(ellipticity) . Thus, the triple of vector fields Vj is a frame. The existence of a 
frame means that the manifold M is parallelizable. 

Conversely, given a frame Vj we uniquely recover the elliptic principal sym- 
bol A\(x, £) via formulae (jlO.ip , (IA.1[) and (|A.2|1 . Thus, a principal symbol is 
equivalent to a frame. 

It is easy to see that the frame elements Vj are orthonormal with respect to 
the metric (|1.32j) . Moreover, the metric can be defined directly from the frame as 

where the repeated frame index j indicates summation over j — 1,2,3. The two 
definitions of the metric, (|1.32[) and (|10.3[) . are equivalent. 

Parallelizability implies orientability. Having chosen a particular orientation 
we define the Hodge star in the standard way. We will use the Hodge star later 
on in this section in order to simplify calculations involving the torsion tensor. 

Note that the topological invariant c introduced in Section Q] in accordance 
with formula (|1.38[) can be equivalently (and more naturally) defined in terms of 
the frame as 

c := sgndetV,". (10.4) 

The crucial new geometric structure is the teleparallel connection. We already 
defined it in Section Q] in accordance with formula (|1.39l) . i.e. via the principal 
symbol. This connection can be equivalently defined via the frame as follows. 
Suppose we have a vector v based at the point y £ M and we want to construct 
a parallel vector u based at the point x G M. We decompose the vector v with 
respect to the frame at the point y. v — c J Vj(y), and reassemble it with the same 
coefficients c J at the point x, defining u :— c J Vj(x). 

We now define the covariant derivative corresponding to the teleparallel con- 
nection. Our teleparallel connection is a special case of an affine connection, so 
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we are looking at a covariant derivative acting on vectors/covectors in the usual 
manner 

= 8v a /dx» + T a ^ v? , V„t»0 = dwp/dx" - w a . (10.5) 

Of course, the above V should not be confused with the V from Section [SJ The 
teleparallel connection coefficients are defined from the conditions 

V M V, a =0, (10.6) 

where the Vj are elements of our frame. Formula (|10.6p gives a system of 27 linear 
algebraic equations for the determination of 27 unknown connection coefficients. 
It is known (see, for example, formula (A2) in [7]), that the unique solution of this 
system is 

r%,3 = v k a (dv kp /dx^, (io.7) 

where 

Vkfj := gfrVk 1 . (10.8) 

The triple of covector fields V^, fc = 1,2, 3, is called the coframe. The frame and 
coframe uniquely determine each other via the relation 

V^Vkv = S jh . (10.9) 

One can check by performing explicit calculations that the teleparallel connec- 
tion has the following two important properties: 

• V Q g ( g 7 = 0, which means that the connection is metric compatible and 

• V Q V/3 — V^V Q = 0, which means that the Riemann curvature tensor is zero. 

The tensor characterising the "strength" of the teleparallel connection is not 
the Riemann curvature tensor but the torsion tensor (|1.40p . The teleparallel con- 
nection is, in a sense, the opposite of the more common Levi-Civita connection: 
the Levi-Civita connection has zero torsion but nonzero curvature, whereas the 
teleparallel connection has nonzero torsion but zero curvature. In our paper we 
distinguish these two affine connections by using different notation for connection 
coefficients: we write the teleparallel connection coefficients as T Q ^ 7 and the Levi- 
Civita connection coefficients (Christoffel symbols) as {^ji see formula (|A.4|) ■ 

Substituting (|10.7[) into (|1.40p we arrive at the following explicit formula for 
the torsion tensor of the teleparallel connection 

T = V j ®dV j , (10.10) 

where the d stands for the exterior derivative. For the sake of clarity we rewrite 
formula (|10.10p in more detailed form, retaining all tensor indices, 

T a M = V^idV^/dx - dV^/dx 1 ) . (10.11) 

As always, the repeated index j appearing in formulae (|10.10p and (jlO.lip indicates 
summation over j = 1,2,3. 
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As pointed out in Section [TJ it is more convenient to work with the rank two 

tensor T defined by formula (|1.41[) rather than with the rank three tensor T. 
Substituting (|10.10l) into (fl~4T"j) we get 

T = Vj®cva\\j , (10.12) 

where 

(curlV^ := (*dV i ) p = \ (dV^e^V^g^. (10.13) 

The remainder of this section is devoted to the proof of formula (|1 .42|) express- 
ing the scalar curvature of the U(l) connection via the torsion of the teleparallel 
connection and the metric. 

We fix an arbitrary point Q £ T'M and prove formula (|1.42p at this point. 
As the LHS and RHS of (|1.42[) are invariant under changes of local coordinates x, 
it is sufficient to prove formula f| 1 .42[) in Riemann normal coordinates, i.e. local 
coordinates such that x = corresponds to the projection of the point Q onto 
the base manifold, g^(0) = 5^ and ^f(O) = 0. Moreover, as the formula we 
are proving involves only first partial derivatives, we may assume, without loss of 
generality, that g^ix) = for all x in some neighbourhood of the origin. Thus, 
it is sufficient to prove formula (|1.42[) for the case of Euclidean metric. 

As both the LHS and RHS of (|1.42[) have the same degree of homogeneity in £, 
namely, — 1, it is sufficient to prove formula (|1 .42|) for £ of norm 1. Moreover, by 
rotating our Cartesian coordinate system we can reduce the case of general £ of 
norm 1 to the case 

£=(0 1). (10.14) 

There is one further simplification that can be made: we claim that it is suffi- 
cient to prove formula (|1.42| for the case when 

VS a (0) = cV, (10.15) 

i.e. for the case when at the point x = the elements of the frame are aligned with 
the coordinate axes; here c = ±1 is the topological invariant defined in accordance 
with formula (|1.38p or, equivalently, in accordance with formula (|10.4II . This claim 
follows from the observation that the LHS of formula (|1.42p is invariant under rigid 
special unitary transformations of the column-function v + (x,£), 

v + i — ^ Rv + , 

where "rigid" refers to the fact that the matrix R £ SU(2) is constant. Of course, 
the column-function Rv + is no longer an eigenvector of the original principal sym- 
bol, but a new principal symbol obtained from the old one by the rigid special 
orthogonal transformation of the frame (|A.27j) with the 3x3 special orthogonal 
matrix O expressed in terms of the 2x2 special unitary matrix R in accordance 
with (|A.28|) . One can always choose the special unitary matrix R so that at the 
point x — the elements of the new frame are aligned with the coordinate axes (in 
fact, there are two possible choices of R which differ by sign). It remains only to 



43 



note the well known fact that the tensor T appearing in the RHS of formula (|1.42l) 
is also invariant under rigid special orthogonal transformations of the frame. 

Having made all the simplifying assumptions listed above, we are now in a 
position to prove formula p. 421) . We give the proof for the case 

c = +l. (10.16) 

There is no need to give a separate proof for the case c = — 1 as the two cases 
reduce to one another by means of the identity (|1.22l) and the observation that 
torsion (|10.10[) is invariant under inversion of the frame. 

Let us calculate the RHS of (|1.42[) first. In view of (j!0.15[) we have, in the 
linear approximation in x, 

f V^[x) Pi 2 (x) V^{x)\ / 1 w 3 (x) -w 2 (x)\ 

V 2 1 {x) V 2 2 (x) V 2 3 (x) = -w 3 {x) 1 w^x) , (10.17) 

^(x) V 3 2 (x) V 3 3 (x)J \w 2 (x) -w\x) 1 / 

where w is some smooth vector-function which vanishes at x = 0. Formula (|10.17jl 
is the standard formula for the linearisation of an orthogonal matrix about the 
identity; see also formula (10.1) in [5]. Note that in Cosserat elasticity literature 
the vector-function w is called the vector of microrotations. Substituting ([10.17D 
into (|10.12p we get, at x = 0, 

T af} = dw /dx a - S aP div to, (10.18) 

which is formula (10.5) from [5|. Here we freely lower and raise tensor indices 
using the fact that the metric is Euclidean (in the Euclidean case it does not 
matter whether a tensor index comes as a subscript or a superscript). Substituting 
(|10.18p and (|10.14[) into the RHS of (fTl^j) we get, at our point Q e T'M, 

1 T a ^ = _ l( dw i/ dx i + Qvf/Q^y ( 10 .i9) 



2 {g^U») m 



Let us now calculate the LHS of (|1.42[) . The equation for the eigenvector 
v ( x >0 of the principal symbol is 



0. (10.20) 



In view of (|10.14p . (|10.15[) and (|10.16[) the (normalised) solution of (|10.20p at our 
point Q e T'M is 

T 

Of course, our v + (x,£) is defined up to the gauge transformation (|1.18l) . (|1.19|) . 
however the LHS of (|1.42p is invariant under this gauge transformation. We now 
perturb equation (|10.20p about the point Q € T'M, that is, about 

x = 0, £ = (0 1) , 
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making use of formula (|10.17|) , which gives us the following equation for the incre- 
ment 8v + of the eigenvector of the principal symbol: 

\ (8v+\ f -w 2 {x) -iw 1 (x)\ fl 

-2J \5v£J \-w 2 (x)+iw 1 (x) J [O 

<5£i - i5&\ fl 



8£i + i6& -26& J V 



= 0, 



or, equivalently, 



8v+ = -(-w 2 (x)+iw 1 (x) + 5£i+i5&)- (10.21) 



Formula (110. 2ip has to be supplemented by the normalisation condition 
||i> + (x, £)|| = 1, which in its linearised form reads 

Re5u+=0. (10.22) 

Formulae (|10.22| and (|10.21j) define Sv + modulo an arbitrary ImSvf, with this 
degree of freedom being associated with the gauge transformation (|1.18|1 . (|1.19|) . 
Without loss of generality we may assume that the gauge is chosen so that 

hn6vf=0. (10.23) 

Combining formulae (|10.22[) . (|10.23[) and (|10.2ip we get 

$ v+ = I ( 2 / \ , ■ 1? N , tt , r t ] • (10.24) 
2 \-W 2 (x) + iw l (x) + 5£i + 18& J 

Recall that the w appearing in this formula is some smooth vector-function which 
vanishes at x = 0. 

Differentiation of (|10.24p gives us 

a " + U ° ' (10.25) 

Formulae (|10.25jl and (|10.26p imply that at our point Q e T'M 

- i{[v+]*, v+} = ~(dw 1 /dx 1 + dw 2 /dx 2 ). (10.27) 

Comparing formulae (|10.19p and (|10.27p and recalling (|10.16p , we arrive at the 
required result (jl.42p . 



dx a 2 \-dw 2 /dx a + idw 1 /dx l 
dv+ 1 /0\ dv+ 1 fo\ dv+ 
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We end this section by writing down an explicit self-contained formula for the 

* 

trace of the tensor T. Note that according to formula (|1.45p , it is only the trace of 

* 

T that we need for our spectral asymptotics. Formulae (|10.12[) and (|10.13[) imply 



trT = Vdet g a ? [V 3l dV ]3 /dx 2 + V i2 dV n /dx 3 + V j3 dV^/dx 1 

- Vji dV j2 /dx 3 - Y j2 dV^/dx 1 - V j3 dV^/dx 2 ] . (10.28) 

Here the coframe Vj is determined from the principal symbol Ai(x,£) in accor- 
dance with formulae (110.11) . (|10.2[) and (|10.8[) or (| 10. 9[) . whereas the metric g is 
determined from the principal symbol Ai(x,^) in accordance with formula (jl.32l) 
or (UTEI]). 

11 Proof of Theorem 11.11 

As Theorem II .11 is an if and only if theorem, our proof comes in two parts. 
11.1 Part 1 of the proof of Theorem 11.11 

Let A be a massless Dirac operator on half-densities. We need to prove that 
a) the subprincipal symbol of this operator, Aj U b(x), is proportional to the identity 
matrix and b) the second asymptotic coefficient of the spectral function, b(x), is 
zero. 

As we have already established the formula for b(x), see (| 1 .45(1 . this part of 
the proof of Theorem 11.11 reduces to proving that the explicit formula for the 
subprincipal symbol of the massless Dirac operator on half-densities is 

A snb (x) = ~(trT(x))l, (11.1) 

where / is the 2x2 identity matrix. 

We give the proof of (|ll.ip for the case (|10. 16|) . There is no need to give a sepa- 
rate proof for the case c = — 1 as the two cases reduce to one another by inversion of 
the frame: the full symbol of the massless Dirac operator on half-densities changes 
sign under inversion of the frame and hence its subprincipal symbol changes sign 
under inversion of the frame, whereas torsion (|10.10[) is invariant under inversion 
of the frame. 

We fix an arbitrary point P € M and prove the identity (jll.ip at this point. 
As the LHS and RHS of (jll.ip are invariant under changes of local coordinates x, 
it is sufficient to check the identity (|ll.ip in Riemann normal coordinates, i.e. lo- 
cal coordinates such that a; = corresponds to the point P, g^iO) = 5^ v and 
^pr(O) = 0. Moreover, as the identity we are proving involves only first partial 
derivatives, we may assume, without loss of generality, that g^vix) — 5^ for all 
x in some neighbourhood of the origin. Furthermore, by rotating our Cartesian 
coordinate system we can achieve (|10.15p . which opens the way to the use, in the 
linear approximation in x, of formula (|10.17l) . 
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Substituting (|10.17j) into (|A~Tj) . we get, in the linear approximation in x, 



,1 



1 + iw 3 



" VI- iw 3 -w 2 ' ai 



1 — iui 1 



C = • 1 2 1 =CT3- (H.2) 

Recall that the w appearing in this formula is some smooth vector-function which 
vanishes at x = 0. 

Substitution of (|11.2[) into (|A.3[) (which coincides with (|A.30p because we as- 
sumed the metric to be Euclidean, g^(x) = 5^) allows us to evaluate the full 
symbol A(x, £) = Ai(x, £) + Aq(x) of the massless Dirac operator on half-densities: 



6s fi - «6 



( w 2 £ 1 -w 1 & iw 3 £i + w 3 & + (-iw 1 - w 2 )£ 3 \ . . 



. i/0 A/0 l\( dw 2 /dx 1 idw 3 /dx 1 \, . . 

Ao(Uj " 4 VI oj 1,1 ^-iS^/ftr 1 -dvP/dx 1 )*-"- [ ' 

Here formula (|11.3p is written in the linear approximation in x, whereas formula 
(| 11 .4[) displays, for the sake of brevity, only one term out of nine (the one corre- 
sponding to a = P = 1 in (|A.3|> ) with the remaining eight terms concealed within 
the dots .... 

Substituting (fTO)) and ([TO]) into (fl~T5]) . we get 

A sub (0) = -i(divw)I. (11.5) 

But, according to (110. 18[) . 

trT(O) = -2divw. (11.6) 
Formulae (fTT3)) , pT5)l and (|10.16p imply formula pTTTjl at x = 0. 



11.2 Part 2 of the proof of Theorem flTTI 

Let A be an operator satisfying assumptions (ll.27[) . (| 1 .3 1 [) and (| 1 . 36[) and such that 
a) the subprincipal symbol of this operator, Aj U b(x), is proportional to the identity 
matrix and b) the second asymptotic coefficient of the spectral function, b(x), is 
zero. We need to prove that A is a massless Dirac operator on half-densities. 

As we have already established the formula for b(x), see (jl.45p . we have, for our 
operator A, the identity (jll.ip . Let Vj be the frame corresponding to the principal 
symbol of the operator A, see formulae (| 10 . 1[) and (|10.2|) . Now, let B be the 
massless Dirac operator on half-densities corresponding to the same frame. Then 
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the principal symbols of the operators A and B coincide. But the subprincipal 
symbols of the operators A and B coincide as well, as in both cases these are 
determined via the frame according to the same formula (for the massless 

Dirac operator B this is the result from subsection 1 A first order differential 
operator is determined by its principal and subprincipal symbols, hence, A = B. □ 

12 Spectral asymmetry 

In this section we deal with the special case when the operator A is differential (as 
opposed to pseudodifferential). No assumptions are made regarding n, m or tr A\. 

Our aim is to examine what happens when we change the sign of the operator. 
In other words, we compare the original operator A with the operator A := —A. 
In theoretical physics the transformation A i— > — A would be interpreted as time 
reversal, see equation (|1.3p . 

It is easy to see that for a differential operator the number to (number of 
equations in our system) has to be even and that the principal symbol has to have 
the same number of positive and negative eigenvalues. In the notation of Section[T] 
this fact can be expressed as to = 2m + = 2mT . 

It is also easy to see that the principal symbols of the two operators, A and A, 
and the eigenvalues and eigenvectors of the principal symbols are related as 





= A 1 (x,-Q, 


(12.1) 






(12.2) 






(12.3) 



whereas the subprincipal symbols are related as 

A suh (x) = -A auh (x). (12.4) 
Formulae JT23||, (QUI , (fl~T7|) . (fl~To) and (Tmj) - (Tl2~4D imply 

a(x) = a(x), b(x) = -b(x). (12.5) 
Substituting (IT23)) into $L8§ and flU]) we get 

a = a, b=-b. (12.6) 

Formulae (|1.7jl and (|12.6|l imply that the spectrum of a generic first order 
differential operator is asymmetric about A = 0. This phenomenon is known in 
differential geometry as spectral asymmetry [TJ [51 [3J 13] ■ 

If we square our operator A and consider the spectral problem A 2 v = A 2 i>, then 
the terms ±b\ n ~ 1 cancel out and the second asymptotic coefficient of the counting 
function (as well as the spectral function) of the operator A 2 turns to zero. This 
is in agreement with the known fact that for an even order semi-bounded matrix 
differential operator acting on a manifold without boundary the second asymptotic 
coefficient of the counting function is zero, see Section 6 of [33] and |24j . 
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The case of the massless Dirac operator is special because, according to Theo- 
rem [TTT1 the spectrum (as well as the spectral function) of this operator is asymp- 
totically symmetric about A = in the two leading terms. However, despite this 
asymptotic symmetry, we believe that for a generic Riemannian 3-manifold the 
spectrum of the massless Dirac operator is asymmetric. In stating this belief we 
are in agreement with the discussion presented on page 1298 of [32) : note that 
in the case of an odd-dimensional manifold the author of [32J refers to the mass- 
less Dirac operator as the Pauli operator. And, of course, our belief that for a 
generic Riemannian 3-manifold the spectrum of the massless Dirac operator is 
asymmetric is closely related to the fact that in dimension 3 the massless Dirac 
operator commutes with the operator of charge conjugation, see formulae (|A.18|) 
and p9) . 

13 Bibliographic review 

To our knowledge, the first publication on the subject of two-term spectral asymp- 
totics for systems was Ivrii's 1980 paper [T3] in Section 2 of which the author stated, 
without proof, a formula for the second asymptotic coefficient of the counting func- 
tion. In a subsequent 1982 paper |15| Ivrii acknowledged that the formula from 
[14] was incorrect and gave a new formula, labelled (0.6), followed by a "proof". In 
his 1984 Springer Lecture Notes [16 Ivrii acknowledged on page 226 that both his 
previous formulae for the second asymptotic coefficient were incorrect and stated, 
without proof, yet another formula. 

Roughly at the same time Rozenblyum [22J also stated a formula for the second 
asymptotic coefficient of the counting function of a first order system. 

The formulae from |TJ] , [T5] and [22] are fundamentally flawed because they are 
proportional to the subprincipal symbol. As our formulae (|1.9p and (|1.24p show, 
the second asymptotic coefficient of the counting function may be nonzero even 
when the subprincipal symbol is zero. This illustrates, yet again, the difference 
between scalar operators and systems. 

The formula on page 226 of [T^ gives an algorithm for the calculation of the 
correction term designed to take account of the effect described in the previous 
paragraph. This algorithm requires the evaluation of a limit of a complicated 
expression involving the integral, over the cotangent bundle, of the trace of the 
symbol of the resolvent of the operator A constructed by means of pseudodiffer- 
ential calculus. This algorithm was revisited in Ivrii's 1998 book, see formulae 
(4.3.39) and (4.2.25) in \T7\. 

The next contributor to the subject was Safarov who, in his 1989 DSc The- 
sis [25], wrote down a formula for the second asymptotic coefficient of the count- 
ing function which was "almost" correct. This formula appears in [35] as for- 
mula (2.4). As explained in Section [TJ Safarov lost only the curvature terms 
I h^{[v Vi]*,v V>}. Safarov's DSc Thesis [25] provides arguments which are 
sufficiently detailed and we were able to identify the precise point (page 163) at 
which the mistake occurred. 

In 1998 Nicoll rederived [20J Safarov's formula (|1.12jl for the principal symbols 
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of the propagator, using a method slightly different from [25], but stopped short 
of calculating the second asymptotic coefficient of the counting function. 

In 2007 Kamotski and Ruzhansky |18| performed an analysis of the propa- 
gator of a first order elliptic system based on the approach of Rozenblyum |22| . 
but stopped short of calculating the second asymptotic coefficient of the counting 
function. 

One of the authors of this paper, Vassiliev, considered systems in Section 6 
of his 1984 paper [33]. However, that paper dealt with systems of a very special 
type: differential (as opposed to pseudodifferential) and of even (as opposed to 
odd) order. In this case the second asymptotic coefficients of the counting function 
and the spectral function vanish, provided the manifold does not have a boundary. 



A The massless Dirac operator 

Let M be a 3-dimensional connected compact oriented manifold equipped with a 
Riemannian metric g a p, ot,/3 = 1,2,3 being the tensor indices. Note that we are 
more prescriptive in this appendix than in the main text of the paper: in the main 
text orientability and existence of a metric emerged as consequences of the way 
we stated the problem, whereas in this appendix they are a priori assumptions. 

We work only in local coordinates with prescribed orientation. 

It is known [311 \W\ that a 3-dimensional oriented manifold is parallelizable, 
i.e. there exist smooth real vector fields Vj, j = 1,2,3, that are linearly indepen- 
dent at every point x of the manifold. (This fact is often referred to as Steenrod's 
theorem.) Each vector Vj(x) has coordinate components Vj a (x), a = 1,2,3. Note 
that we use the Latin letter j for enumerating the vector fields (this is an anholo- 
nomic or frame index) and the Greek letter a for enumerating their components 
(this is a holonomic or tensor index). The triple of linearly independent vector 
fields Vj, j — 1, 2, 3, is called a frame. Without loss of generality we assume further 
on that the vector fields Vj are orthonormal with respect to our metric: this can 
always be achieved by means of the Gram-Schmidt process. 

Define Pauli matrices 

a a (x) :=s j Vj a (x), (A.l) 

where 

s 1 := (J £)=-!, , 2 := (j "») = s 2 , s 3 := (j = s 3 . (A.2) 

In formula (jA.lj) summation is carried out over the repeated frame index j = 1, 2, 3, 
and a = 1, 2, 3 is the free tensor index. 

The massless Dirac operator is the matrix operator 

— (£+M£ + {£H)- <"> 

where summation is carried out over a, ft, 7 = 1,2,3, and 

/ £ 1 ._ 1 n 0S f 9 9-yS dg aS dg a ^\ 

IcryJ " 2 9 \dx a + dx-r dx 5 J [ ' 
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are the Christoffel symbols. Here and throughout this appendix we raise and lower 
tensor indices using the metric. Note that we chose the letter U W for denoting 
the massless Dirac operator because in theoretical physics literature it is often 
referred to as the Weyl operator. 

Formula (|A.3|) is the formula from |10j , only written in matrix notation (i.e. with- 
out spinor indices). Note that in the process of transcribing formulae from [TP] 
into matrix notation we used the identity 

ea a e = (a a ) T , (A.5) 

a = 1, 2, 3, where 

«=-(!o) <"> 

is the 'metric spinor'. The identity (|A.5|) gives a simple way of raising/lowering 
spinor indices in Pauli matrices in the non-relativistic (a ^ 0) setting. 

Physically, our massless Dirac operator (|A.3j) describes a single neutrino living 
in a 3-dimensional compact universe M. The eigenvalues of the massless Dirac 
operator are the energy levels. 

Observe that the sign of dct Vj a is preserved throughout the connected oriented 
manifold M. Having detVj" > means that our frame has positive orientation 
(relative to the prescribed orientation of local coordinates) and det Vj a < means 
that our frame has negative orientation. Accordingly, we say that our massless 
Dirac operator (|A.3|) has positive/negative orientation depending on the sign of 
det Vj a . Of course, the transformation W H > — W changes the orientation of the 
massless Dirac operator. 

The massless Dirac operator (|A.3|) acts on columns v — (vi v%) of complex- 
valued scalar functions. In differential geometry this object is referred to as a 
(Weyl) spinor so as to emphasise the fact that v transforms in a particular way 
under transformations of the orthonormal frame V. However, as in our exposition 
the frame V is assumed to be chosen a priori, we can treat the components of the 
spinor as scalars. This issue will be revisited below when we state Property 4 of 
the massless Dirac operator. 

We now list the main properties of the massless Dirac operator. 



Property 1. The massless Dirac operator is invariant under changes of local 
coordinates x, i.e. it maps 2-columns of smooth scalar functions M —> C 2 to 2- 
columns of smooth scalar functions M — > C 2 regardless of the choice of local 
coordinates. 

In order to establish this property we examine separately the two operators 

(A.7) 



and 



dx a 



^I^ + Pk] (A. 



dx a \ a^j 



appearing in formula (| A.3|) . 
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Let us act with the differential operator (|A.7[) on a 2-column u : M — > C 2 
of smooth scalar functions. Then is a column-valued covector (i.e. pair of 
gradients), a a is a matrix- valued vector, so matrix multiplication combined with 
contraction in a gives a column- valued scalar. Thus, the operator (|A.7[) is invariant 
under changes of local coordinates. 

As to the multiplication operator (|A.8|) . its invariance follows from the obser- 
vation that ^ff^ { A } <j7 ) * S a ma t r i x-varue d tensor. 

Property 2. The massless Dirac operator is formally self-adjoint (symmetric) 
with respect to the inner product 



v*w \J det g a p dx (A. 9) 

M 



on 2-columns of smooth scalar functions v,w : M 
Indeed, the adjoint operator is 



(A.IO) 



Comparing formulae (jA.3[) and (jA.lOp we see that in order to prove formal self- 
adjointness we need to show that 



£+{£W"- + -"(£ + {£H 



d 



%/det g KX \dx 



— ^tetg^a"). (A.ll) 



We fix an arbitrary point P € M and prove the identity (|A.11[) at this point. 
In view of Property 1, it is sufficient to check the identity (jA.llj) in Riemann 
normal coordinates, i.e. local coordinates such that x = corresponds to the 
point P, g M! ,(0) = 5^ and ^pr(O) = 0. Moreover, as the identity we are proving 
involves only first partial derivatives, we may assume, without loss of generality, 
that g^vix) — S^v for all x in some neighbourhood of the origin. Thus, the problem 
has been reduced to proving that variable (i.e. dependent on a;) Pauli matrices in 
Euclidean space satisfy the identity 

(£)•"■♦-■'(£)-<(£)■ 

Note that in (|A.12p we made all the tensor indices upper, using the fact that the 
metric is Euclidean (in the Euclidean case it does not matter whether a tensor index 
comes as a subscript or a superscript). Of course, we still retain the convention of 
summation over repeated indices. 

In order to prove (|A.12[) we recall the basic identity for Pauli matrices which 
in the Euclidean case reads 

irV + a" a* = 2IS" U , (A.13) 
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where / is the 2x2 identity matrix. (For a general metric one would have written 
the above formula with g^ u instead of S^ u .) Formula (|A.13[) implies 

o-^ = 3/, (A.14) 

er'W 1 = -cr K , (A.15) 
6(0*0" + a u a fl )/dx x = 0. (A.16) 
Using formulae (|A.13|) — (|A.16|) we get 

d<j«\ , fda a \ _ B _ B , _ B fdo<*\ _ 3 , _ g fda 



= a?^ ^— + ^- <^ ^ 



da;"/ \dx a J \dx a J \dx a 



\dx a J \dx a 



a" 



da a \ fda a \ ( da a \ ( da a 



dx a J \dx a J \dx a J \dx a 

\dx a J \dx a J \dx a 

Comparing the left- and right-hand sides of (|A.17[) we arrive at (I A. 12|) . 



(A.17) 



Property 3. The massless Dirac operator W commutes 

C(Wv) = WC(v) (A.18) 

with the antilinear map 

v^C(v):=ev. (A.19) 

Here the map (IA. 19|) acts on columns v = (yi V2) T of complex- valued scalar 
functions, with e being the 'metric spinor' defined in accordance with (|A.6j) . The 
commutativity property (|A. 18|) follows from the explicit formula for the massless 
Dirac operator (|A.3|) and the identity ea a = —a a e, a = 1,2,3, the latter being a 
consequence of formula (|A.5I) . 

Formula (|A.18[) implies that v is an eigenfunction of the massless Dirac oper- 
ator corresponding to an eigenvalue A if and only if C(v) is an eigenfunction of 
the massless Dirac operator corresponding to the same eigenvalue A. Hence, all 
eigenvalues of the massless Dirac operator have even multiplicity. Moreover, any 
eigenfunction v and its 'partner' C(i>) make the same contribution to the spectral 
function (|1.4[) at every point x of the manifold M. 

We do not use the commutativity property (|A.18p of the massless Dirac oper- 
ator in the current paper. 
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The antilinear operator (|A.19[) is, of course, the charge conjugation operator 
which we already encountered in Section [TJ see formula f| 1 . 28[) . The difference 
between the arguments presented in this appendix and those in Section Q] is that 
in this appendix we deal with the differential operator, whereas in Section Q] we 
dealt with the principal symbol. This leads to opposite commutation proper- 
ties: the charge conjugation operator commutes with the Weyl operator but it 
anticommutes with its principal symbol. The source of this difference is the i 
appearing in the RHS of formula (|A.3I) . 

Property 4. This property has to do with a particular behaviour under SU(2) 
transformations. Let R : M — > SU(2) be an arbitrary smooth special unitary 
matrix-function. Let us introduce new Pauli matrices 

a a := Ra a R* (A.20) 

and a new operator W obtained by replacing the a in (|A.3I) by a . It turns out 
(and this is Property 4) that the two operators, W and W, are related in exactly 
the same way as the Pauli matrices, a and a, that is, 

W = RWR*. (A.21) 

In order to prove formula (|A.2ip we write down the operator W explicitly and 
rearrange terms: 

W := -iRa a R* ( * + W ( B ^f^ + ( ^ R<r~< R* 
\dx a 4 ' \ dx a [asy 

= -iRa a -r-R* + iRa a - 



i 



dx a dx a 



. Rv a aR l— + i f \a< )R* - -Ra a (jRR* t^^R* + Ra 
4 p \dx a \ajj J 4 p \dx a dx a 

= RWR* + iRa a ^ - -Ra a asR* ( l^iT + R<t? — 



dx a 4 p \dx a dx a 
Hence, proving (|A.2I[) reduces to proving that 

In order to prove formula (|A.22ll it is sufficient to show that 

n * dR 8 n * g dR * * dR * 
(JrR - — cr' j R* + apaP- — = 4- — 

' dx a 1 dx a dx a 

which, in turn, in view of the identity a per 13 — 31 (we already used it in the special 
case of Euclidean metric, see formula (jA.f 4[) ). is equivalent to proving that 

af}R *™P = ^L R . (A.23) 

' dx a dx a 
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The fact that the matrix function R is special unitary implies that at every point 
x of the manifold M and for every index a — 1, 2, 3 the matrix R*-§§s is trace-free 
anti-Hcrmitian, which, in view of the identity apa 1 a? = —a 1 (we already used it 
in the special case of Euclidean metric, see formula (|A.15[) ). implies that formula 
(1A.23[) can be equivalently rewritten as 

-IT™™R. (A.24) 

dx a dx a 

But formula (|A.24[) is an immediate consequence of the identity R*R = I. 

Having proved Property 4, let us examine the geometric meaning of the trans- 
formation (|A.20p . Let us expand the new Pauli matrices a with respect to the 
basis (|A.2|) : 

a a (x) = s j Vj a (x). (A.25) 

Formulae (|A.1[) . (|A.25|) and (|A.20[) give us the following identity relating the new 
vector fields V^ and the old vector fields V 3 : 

Rs k R*V k = s j Vj . (A.26) 

Resolving (|A.26[) for Vj we get 

Vj = Oj k V k , (A.27) 
where the real scalars Oj k are given by the formula 

Oj k = ^tT( Sj Rs k R*). (A.28) 

Note that in writing formulae (|A.26|I and (|A.27|) we chose to hide the tensor index, 
i.e. we chose to hide the coordinate components of our vector fields. Say, formula 
(|A.27[) written in more detailed form reads Vj a = Oj k Vk a ■ 

The scalars (|A.27[) can be viewed as elements of a real 3x3 matrix-function O 
with the first index, j, enumerating rows and the second, k, enumerating columns. 
It is easy to check that this matrix-function O is special orthogonal. Hence, the 
new vector fields Vj are orthonormal and have the same orientation as the old 
vector fields Vj . We have shown that the transformation (|A.20[) has the geometric 
meaning of switching from our original oriented orthonormal frame Vj to a new 
oriented orthonormal frame Vj . 

Formula (|A~28)) means that the special unitary matrix R is, effectively, a square 
root of the special orthogonal matrix O. It is easy to see that for a given matrix 
O G SO(3) formula (|A.28|) defines the matrix R £ SU(2) uniquely up to sign. This 
observation allows us to view the issue of the geometric meaning of the transfor- 
mation (|A.20[) the other way round: given a pair of orthonormal frames, Vj and Vj, 
with the same orientation, we can recover the special orthogonal matrix- function 
0(x) from formula (|A.27|) and then attempt finding a smooth special unitary 
matrix-function R(x) satisfying (|A.28[) . Unfortunately, this may not always be 
possible due to topological obstructions. We can only guarantee the absence of 
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topological obstructions when the two frames, Vj and Vj, are sufficiently close to 
each other, which is equivalent to saying that we can only guarantee the absence 
of topological obstructions when the special orthogonal matrix- function 0[x) is 
sufficiently close to the identity matrix for all x G M. 

We illustrate the possibility of a topological obstruction by means of an explicit 
example. 

Example A.l Consider the unit torus T 3 parameterized by cyclic coordinates x a , 
a = 1,2,3, of period 2tt. The metric is assumed to be Euclidean. Define the 
orthonormal frame as 

^cos/c 3 a; 3 \ /— sinfc 3 a; 3 

\\" | sin/,;;.r :; I . IV' ( ms/.-;;.r ;i ]. IV | | . (A. 2!) ) 

where fc 3 G Z is a parameter. Let W be the massless Dirac operator corresponding 
to the frame \A . 29\) with some even fc 3 and let W be the massless Dirac operator 
corresponding to the frame HA . 29\) with some odd fc 3 . We claim that there does 
not exist a smooth matrix-function R : T 3 — > SU(2) which would give HA.28\) . 
where 0(x) is the special orthogonal matrix-function defined by formula &A.27\) . 
To prove this, it is sufficient to show that the two operators, W and W , have 
different spectra. Straightforward separation of variables shows that any half-even 
integer (positive or negative) is an eigenvalue of W but is not an eigenvalue of 
W . What happens in this example is that a special unitary matrix-function R(x) 
satisfying \A . 28\) can be defined locally but not globally: if we try to construct 
R(x 3 ) moving along the circumference of the torus x 3 € (— 7T, 71") we end up with a 
discontinuity, lim R(x 3 ) = — lim R(x 3 ). 

x 3 — 7r+ a; 3 — >tt~ 

In fact, one can generalise Example IA.1I by introducing rotations in three dif- 
ferent directions, which leads to eight genuinely distinct parallelizations. See also 
[25] page 524. 

Let us emphasise that the topological obstructions we were discussing have 
nothing to do with Stiefel-Whitney classes. We are working on a parallelizable 
manifold and the Stiefel-Whitney class of such a manifold is trivial. The topologi- 
cal issue at hand is that our parallelizable manifold may be equipped with different 
spin structures. 

We say that two massless Dirac operators, W and W, are equivalent if there 
exists a smooth matrix-function R : M — > SU(2) such that the corresponding Pauli 
matrices, a a and a a , are related in accordance with (IA.20[) . In view of Property 4 
(see formula (IA.21[) ') all massless Dirac operators from the same equivalence class 
generate the same spectral function (jl.4l) and the same counting function (11.51) , so 
for the purposes of our paper viewing such operators as equivalent is most natural. 

As explained above, there may be many distinct equivalence classes of massless 
Dirac operators, the difference between which is topological. Studying the spectral 
theoretic implications of these topological differences is beyond the scope of our 
paper. The two-term asymptotics (jl.6|) and (jl.7|) derived in the main text of our 
paper do not feel this topology. 
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In theoretical physics the SU(2) freedom involved in defining the massless Dirac 
operator is interpreted as a gauge degree of freedom. We do not adopt this point 
of view (at least explicitly) in order to fit the massless Dirac operator into the 
standard spectral theoretic framework. 

We defined the massless Dirac operator (IA.3|) as an operator acting on 2- 
columns of scalar functions, i.e. on 2-columns of quantities which do not change 
under changes of local coordinates. This necessitated the introduction of the den- 
sity i^/det g a p in the formula (|A.9[) for the inner product. In spectral theory it is 
more common to work with half-densities. Hence, we introduce the operator 

W 1/2 := (det 5KA ) 1 / 4 I¥(det 5 ^)- 1/4 (A.30) 

which maps half-densities to half-densities. We call the operator (|A.30[) the mass- 
less Dirac operator on half- densities. 
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